MULTIPLE HAMILTONIAN STRUCTURE 
OF BOGOYAVLENSKY-TODA LATTICES 



Pantelis A. Damianou 



Department of Mathematics and Statistics 
University of Cyprus 
P. O. Box 20537, 1678 Nicosia, Cyprus 
Email: damianou@ucy.ac.cy 



ABSTRACT 

This paper is mainly a review of the multi-Hamiltonian nature of Toda and generalized Toda 
lattices corresponding to the classical simple Lie groups but it includes also some new results. 
The areas investigated include master symmetries, recursion operators, higher Poisson brack- 
ets, invariants and group symmetries for the systems. In addition to the positive hierarchy 
we also consider the negative hierarchy which is crucial in establishing the bi-Hamiltonian 
structure for each particular simple Lie group. Finally, we include some results on point and 
Noether symmetries and an interesting connection with the exponents of simple Lie groups. 
The case of exceptional simple Lie groups is still an open problem. 
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1 INTRODUCTION 



In this paper we review the bi-Hamiltonian and multiple Hamiltonian nature of the Toda 
lattices corresponding to simple Lie groups. These are systems that generalize the usual 
finite, non-periodic Toda lattice (which corresponds to a root system of type A n ). This 
generalization is due to Bogoyavlensky [1]. These systems were studied extensively in [2] 
where the solution of the system was connected intimately with the representation theory of 
simple Lie groups. There are also studies by Olshanetsky and Perelomov [3] and Adler, van 
Moerbeke [4]. We will call such systems the Bogoyavlensky- Toda lattices. 
We begin with the following more general definition which involves systems with exponential 
interaction: Consider a Hamiltonian of the form 

i m 

# = ^(p,p) + E e(Vl,q) > (!) 

z i=l 

where q = (qi, . . . , qx), p = (pi, • • • ,Pn), Vi, . . . , v m are vectors in and ( , ) is the 
standard inner product in R^. The set of vectors A = {vi, . . . , v m } is called the spectrum 
of the system. 

In this paper we limit our attention to the case where the spectrum is a system of simple 
roots for a simple Lie algebra Q. In this case m = I = rankC?. It is worth mentioning that 
the case where m, N are arbitrary is an open and unexplored area of research. The main 
exception is the work of Kozlov and Treshchev [5] where a classification of system (1) is 
performed under the assumption that the system possesses N polynomial (in the momenta) 
integrals. We also note the papers by Ranada [6] and Annamalai, Tamizhmani [7]. Such 
systems are called Birkhoff integrable. For each Hamiltonian in (1) we associate a Dynkin 
type diagram as follows: It is a graph whose vertices correspond to the elements of A. Each 
pair of vertices v iy Vj are connected by 

{Vi,Vi)(Vj,Vj) 

edges. 

Example: The usual Toda lattice corresponds to a Lie algebra of type A^-i- In other 
words m — I — N — 1 and we choose A to be the set: 

v 1 = (1,-1,0,..., 0), ,v N -i = (0,0,..., 0,1,-1) . 

The graph is the usual Dynkin diagram of a Lie algebra of type Atv-i- The Hamiltonian 

becomes 

N i N-l 

H( Ql , ...,q N , Pl ,... ,p N ) = £ -p? + £ e*-*" , (2) 

i=i z i=i 

which is the well-known classical, non-periodic Toda lattice. 

It is more convenient to work, instead in the space of the natural (q,p) variables, with the 
Flaschka variables (a, b) which are defined by: 



i — 1, 2, . . . , m 
* = l,2,...,iV 



(3) 



We end-up with a new set of polynomial equations in the variables (a, b). One can write 
the equations in Lax pair form (at least this is well-known for spectrum corresponding to 
simple Lie algebras); see for example [8]. The Lax pair (L(t), B(t)) in Q can be described in 
terms of the root system as follows: 

i i 

L(t) = k{t)h ai + <k(t)(xai + z-cj , 
i=i i=i 

i 

B{t) =^2a i (t)(x a . -x_ Q J . 
i=i 

As usual h ai is an element of a fixed Cartan subalgebra and x ai is a root vector corresponding 
to the simple root aj. The Chevalley invariants of Q provide for the constants of motion. 
We will describe them separately for each case. 

In this paper we begin with a review of the An Toda system. The theory and bi-Hamiltonian 
structure for this case is well-developed and the results are well-known. We present a review 
of the results in sections 3-5 and then, for the remaining part of the paper, we deal exclusively 
with the other classical simple Lie algebras of type B^, and D^. We will demonstrate 
that these systems are bi-Hamiltonian and then illustrate with some small dimensional 
examples, namely B 2 , C 3 and D 4 . 

The multi-Hamiltonian structure of the Toda lattice was established in [9]. For the remaining 
Bogoyavlensky-Toda systems it was established recently in several papers. The results for 
Bn Bogoyavlensky-Toda lattice were computed in [10] in Flaschka coordinates, and in [11] 
in (q,p) coordinates. The Cn case is in [12] in Flaschka coordinates, and [11] in natural 
(q,p) coordinates. The D N case was settled in [13]. The bi-Hamiltonian structure of these 
systems was established recently in [14]. The negative Toda hierarchy was constructed in 
[15] and it was crucial in establishing the bi-Hamiltonian formulation in [14]. 
The construction of the bi-Hamiltonian pair may be summarized as follows: 
Define a recursion operator 1Z in (a, b) space by finding a second bracket, 7r 3 , and inverting 
the initial Poisson bracket tx\. Define the negative recursion operator Af by inverting the 
second Poisson bracket ty^. This recursion operator is the inverse of the operator TZ. Finally, 
define a new rational bracket 7r_i by 7r_i = A/Vi = ttiTT^tti. We obtain a bi-Hamiltonian 
formulation of the system: 

where Hi = jtrL*. The brackets 7Ti and 7r_i are compatible and Poisson. 
There is also an interesting connection with the exponents of the corresponding Lie group. 
For example, in the case of Dn there is a sequence of invariants H2, H4, . . . , of even degree 
and an additional invariant of degree N. Let Xi denote the Hamiltonian vector field generated 
by Hi and let Zq denote a conformal symmetry. Then we have 

[Zo,Xj] = f(j)Xj ■ 

The values of f(j) corresponding to independent Xj generate all the exponents except one. 
When Zq acts on the Hamiltonian vector field XPi where P is the invariant corresponding 
to the Pfafnan of the Jacobi matrix, we obtain the last exponent N — 1. For example, in the 
case of D 5 the exponents are 1, 3, 5, 7, and 4 . The independent invariants are if 2, H4, H 6 , 
H$ and P 5 where H 2 i = h Tr L 2t and P 5 = Vdet L. We obtain 



[Z ,X2] 

[Zo,Xe] 
[Zo,Xs] 
[^o,Xp 5 ] 

In other words, the coefficients on the right 1 
Lie group of type D 5 . 



= X2 

= 3x4 
= 5X6 
= 7xs 
= 4 X p 5 • 

nd side are precisely the exponents of a simple 



2 BACKGROUND 

In this section we review the necessary background from Poisson and symplectic geometry, 
bi-Hamiltonian systems, master symmetries and recursion operators. 



2.1 Schouten bracket 

We list some properties of the Schouten bracket following [16], [17], [18]. Let M be a C°° 
manifold, N = C°°(M) the algebra of C°° real-valued functions on M. A contravariant, 
antisymmetric tensor of order p will be called a p-tensor for short. These tensors form a 
superspace endowed with a Lie-superalgebra structure via the Schouten bracket. 

The Schouten bracket assigns to each p-tensor A, and g-tensor B, a {p+q — l)-tensor, denoted 
by [A, B]. For p = 1 we have [A, B] = L^B where La is the Lie-derivative in the direction 
of the vector field A. The bracket satisfies 



[A,B] = (-ir[B,A] . (4) 

ii) If C is a r-tensor 

C],A] + (-ir[[C, A],B] + (-ir[[A, B],C] = 0. (5) 

Hi) 

[A, B A C] = [A, B] A C + {-l) pq+q B A [A, C] . (6) 



2.2 Poisson Manifolds 

We review the basic definitions and properties of Poisson manifolds following [16], [18], [19]. 
A Poisson structure on M is a bilinear form, called the Poisson bracket { , } : N x N — > N 
such that 

{f,g} = -{g,f} (7) 

ii) 

{/, {<?, h}} + {g, {h, /}} + {h, {/, = (8) 
= + (9) 



Properties i) and ii) define a Lie algebra structure on N. ii) is called the Jacobi identity 
and Hi) is the analogue of Leibniz rule from calculus. A Poisson manifold is a manifold M 
together with a Poisson bracket { , }. 

To a Poisson bracket one can associate a 2-tensor % such that 

{/ 5 9} — ( 7r ) df A dg) . (10) 

Jacobi's identity is equivalent to the condition [n, n] = where [ , ] is the Schouten bracket. 
Therefore, one could define a Poisson manifold by specifying a pair (M, ir) where M is a 
manifold and n a 2-tensor satisfying [n, n] — 0. In local coordinates (xi,x 2 , ■ ■ ■ ,x n ), n is 
given by 

and 

{/,«,} = (7r,d/ A ^> = E^A||. (12) 

In particular {xi,Xj} = 7Tij(x). Knowledge of the Poisson matrix (ir^) is sufficient to define 
the bracket of arbitrary functions. The rank of the matrix (71^) at a point x G M is called 
the rank of the Poisson structure at x. 

A function F : M 1 — > M 2 between two Poisson manifolds is called a Poisson mapping if 

{foF,goF} 1 = {f,g} 2 oF (13) 

for all /, g e C°°(M2). In terms of tensors, F*7Ti = 7r 2 . Two Poisson manifolds are called 
isomorphic, if there exists a diffeomorphism between them which is a Poisson mapping. 

The Poisson bracket allows one to associate a vector field to each element / e N. The vector 
field Xf is defined by the formula 

X/(<7) = {/,<?}■ (14) 
It is called the Hamiltonian vector field generated by /. In terms of the Schouten bracket 

Xf = WJ]. (15) 

Hamiltonian vector fields are infinitesimal automorphisms of the Poisson structure. These 
are vector fields X satisfying L x ^ = 0. In the case of Hamiltonian vector fields we have 

L X} k = [tt.x/] = [f, [*,/]] = -2[[7r,7r],/] = . (16) 

The Hamiltonian vector fields form a Lie algebra and in fact 

[Xf,Xg] = X{f,g} ■ (17) 

So, the map / — > is a Lie algebra homomorphism. 
The Poisson structure defines a bundle map 

7T* : T*M -> TM (18) 

such that 

**m=xf- (19) 



The rank of the Poisson structure at a point x G M is the rank of n* : T*M — > T X M. 
Throughout this paper we use the symbol tt to denote a Poisson tensor but occasionally we 
use the same symbol to denote the matrix of the components of the tensor (i.e. the Poisson 
matrix). The same convention applies for the recursion operators. 

The functions in the center of N are called Casimirs. It is the set of functions / so that 
= f° r an 9 £ N. These are functions which are constant along the orbits of 
Hamiltonian vector fields. The differentials of these functions are in the kernel of n*. In 
terms of the Schouten bracket a Casimir satisfies [n, f] = 0. 

Given a function /, there is a reasonable algorithm for constructing a Poisson bracket in 
which / is a Casimir. One finds two vector fields X 1 and X 2 such that L Xl f = L X2 f = 0. 
If in addition X±, X 2 and [X 1? X 2 ] are linearly dependent, then X\ A X 2 is a Poisson tensor 
and / is a Casimir in this bracket. In fact 

[/, X 1 A X 2 ] = [f, X,} A X 2 - X 1 A [/, X 2 ]=0. (20) 

More generally, there is a formula due to Flaschka and Ratiu which gives locally a Pois- 
son bracket when the number of Casimirs is 2 less than the dimension of the space. Let 
/x, / 2 , . . . , f r be functions on R r+2 . 
Then the formula 

u>{g, h} — dfx A . . . A df r A dg A dh (21) 

where u is a non- vanishing r + 2 form, defines a Poisson bracket on R r+2 and the functions 
/x, . . . , f r are Casimirs. For more details on this formula, see [20]. 

Multiplication of a Poisson bracket by a Casimir gives another Poisson bracket. Suppose 
[tt, tt] = and [tt, /] = 0. Then 

[fir, fn] = f A [/, tt] A tt + / A tt A [tt, /] + f 2 [n, tt] = . (22) 



2.3 Symplectic and Lie Poisson manifolds 

The most basic examples of Poisson brackets are the symplectic and Lie-Poisson brackets. 

i) Symplectic manifolds: A symplectic manifold is a pair (M 2n , oo) where M 2n is an even 
dimensional manifold and uj is a closed, non-degenerate 2-form. The associated isomorphism 

fi-.TM -> T*M (23) 

extends naturally to a tensor bundle isomorphism still denoted by /i. Let A = /i" 1 , f E N 
and let x/ — be the corresponding Hamiltonian vector field. The symplectic bracket 

is given by 

{f,g} = u>(Xf,Xg) ■ (24) 

In the case of R 2n , according to a Theorem of Darboux, there are coordinates (xi, . . . , x n , yi, . . . , y n ), 
so that 

n 

uj = dxi A dyi (25) 

i=l 



and the Poisson bracket is the standard symplectic bracket on R 2n . 

ii) Lie-Poisson : Let M = Q* where Q is a Lie algebra. For aG§, define the function $ a 
on Q* by 

$M = (a,ti) (26) 

where \i G Q* and ( , ) is the pairing 
between Q and Q*. Define a bracket on Q* by 

{$ a ,$ 6 } = $ M . (27) 

This bracket is easily extended to arbitrary C°° functions on Q*. The bracket of linear func- 
tions is linear and every linear bracket is of this form, i.e., it is associated with a Lie algebra. 
Therefore, the classification of linear Poisson brackets is equivalent to the classification of 
Lie algebras. 



2.4 Local theory 

In his paper [19] A. Weinstein proved the so-called "splitting theorem" , which describes the 
local behavior of Poisson manifolds. 

Theorem 1 Let x be a point in a Poisson manifold M. Then 

near Xq, M is isomorphic to a product SxN where S is symplectic, N is a Poisson manifold, 
and the rank of N at x is zero. 

S is called the symplectic leaf through x$ and N is called the transverse Poisson structure 
at xq. N is unique up to isomorphism. So, through each point x$ passes a symplectic leaf 
S XQ whose dimension equals the rank of the Poisson structure on M aX xq. The bracket on 
the transverse manifold N xo can be calculated using Dirac's constraint bracket formula. 

Theorem 2 Let x be a point in a Poisson manifold M 

and let U be a neighborhood of Xq which is isomorphic to a product 

S x N as in Weinstein 's splitting Theorem. Let p i} i — 1, . . . , 2n be functions on U such that 

N = {x eU\pi(x) = constant} . (28) 

Denote by P = Pij = {pi,Pj} and by P u 

the inverse matrix of P. Then the bracket formula for the 

transverse Poisson structure on N is given as follows: 

2n 

{F,G} N (x) = {F,G} M (x) + ^2{F,p l } M (x)P l ^x){G,p,} M (x) (29) 

i,3 

for all x G N, where F , G are functions on N and F , G 

are extensions of F and G to a neighborhood of M. Dirac's formula depends 

only on F, G, but not on the extensions F, G. 



When 11 is an element of Q*, where Q is a semi-simple Lie algebra, Cushman and Roberts 
proved that, in suitable coordinates, the transverse structure is polynomial; see [21] and [22]. 



2.5 Cohomology 

Cohomology of Lie algebras was introduced by Chevalley and Eilenberg in [23]. Let Q be a 
Lie algebra and let p be a representation of Q with representation space V. A (/-linear skew- 
symmetric mapping of Q into V will be called a g-dimensional IZ-cochain. The g-cochains 
form a space C q (G, V). By definition, C°(G, V) = V. 

We define a coboundary operator 5 = 5 q : C q (Q, V) — > C" ?+1 (^, V) by the formula 
(5/)(x ,...,Xg) = J2i=o(-l) q p(xi)f(xo,---,Xi,...,x q ) + 

(30) 

Si<j ( 1) "'/( [•£? ) -^j] j *^0i ■ ■ ■ , . . . , Xj , . . . , Xq) , 

where / G C q (Q,V) and x , . . . , x g G C?. As can be easily checked S q+1 o S q = so that 
{C q (Q, V), 5 q } is an algebraic complex. Define Z q (Q, V) the space of g-cocycles as the kernel 
of S : C q — > C" ?+1 and the space -B 9 ({?, V) of g-coboundaries as the image 5C 9-1 . Since M = 
we can define 

"<*v>-*g#. (3D 

Lichnerowicz [16] considers the following cohomology defined on the tensors of a Poisson 
manifold. Let (M, tc) be a Poisson manifold. If we set B = C = n in (5) we get 

[n,[n,A]] = (32) 

for every tensor A. Define a coboundary operator d n which assigns to each p-tensor A, a 
(p + l)-tensor d n A given by 

d n A = -[n, A] . (33) 

We have d^A = [ty, [tt, A]} = and therefore d n defines a cohomology. An element A is a 
p-cocycle if [tt, A] = 0. An element B is a p-coboundary if B = [ty, C], for some (p — l)-tensor 
C. Let 

Z n (M,n) = {AeT n \ [tt,A} = 0} (34) 

and 

B n (M,n) = {B\B=[7i,C], CGT n _i}. (35) 

The quotient 

v ' ; B n (M, ty) v ; 

is the nth cohomology group. 

Let £ be a Lie algebra and consider the Lie- Poisson manifold Q*. Define a representation p 
of g with values in C°°{g*) by 

p(*o/ = E4-|r (37) 



where Xi denotes coordinates on Q* and at the same time elements of a basis for Q. In other 
words, p{xj)f = {xi, /}, where the bracket is the Lie-Poisson bracket on Q* . We denote the 
nth cohomology group of Q with respect to this representation by 

H n (g,c°°(g*)) . (38) 

We have the following result: 
Theorem 3 

H n (g*,7v) = w\g, c°°(g*)) . (39) 

The proof can be found in [24] or [10]. 
2.6 Bi Hamiltonian systems 



Proposition 1 Let (M,ivi), (M, 7r 2 ) two Poisson structures on M. The following are equiv- 
alent: 

i) 7Ti + 7r 2 is Poisson. 

ii) [ni,7Y 2 ] = 0. 

Hi) d 7ri d W2 d W2 d 7ri . 

IV) TTi G Z 2 (M,7T 2 ), 7T 2 G Z 2 {M,Tl l ). 

Two tensors which satisfy the equivalent conditions are said to form a Poisson pair on M. 
The corresponding Poisson brackets are called compatible. 

Lemma 1 Suppose ~K\ is Poisson and n 2 = L x ^i = —d ni X for some vector field X. Then 
iii is compatible with tt 2 . 

Proof: 

[7Ti,7r 2 ] = [tti, -[n^X]} = -d Wl d Wl X = . 

■ 

If 7Ti is symplectic, we call the Poisson pair non-degenerate. If we assume a non-degenerate 
pair we make the following definition: The recursion operator associated with a non-degenerate 
pair is the (1, l)-tensor 1Z defined by 

n = n 2 7i^\ (40) 

A bi-Hamiltonian system is defined by specifying two Hamiltonian functions Hi, H 2 satisfy- 
ing 

X = TTi VH 2 = 7T 2 VHi . (41) 

We have the following result due to Magri [25] : 



Theorem 4 Suppose that we have a non-degenerate bi-Hamiltonian system on a manifold 
M , whose first cohomology group is trivial. Then, there exists a hierarchy of mutually com- 
muting functions H 1 , H 2 , . . ., all in involution with respect to both brackets. If we denote by 
Xi the Hamiltonian vector field generated by Hi with respect to the initial bracket 7Ti then 
the Xi generate mutually commuting bi-Hamiltonian flows, satisfying the Lenard recursion 
relations 

Xi+j = m VH, , (42) 
where 7Tj = TZ 1 ~ 1 tc 1 are the higher order Poisson tensors. 

For further information on bi-Hamiltonian systems relevant to Toda type systems see [26], 
[27], [28], [29]. 



2.7 Master Symmetries 

We recall the definition and basic properties of master symmetries following Fuchssteiner 
[30]. Consider a differential equation on a manifold M defined by a vector field x- We are 
mostly interested in the case where x is a Hamiltonian vector field. A vector field Z is a 
symmetry of the equation if 

[Z,X} = 0. (43) 
If Z is time dependent, then a more general condition is 

^ + [Z,x}=0. (44) 
A vector field Z is called a master symmetry if 

[[Z, X ],X}=0, (45) 

but 

[Z, X ]^0. (46) 

Master symmetries were first introduced by Fokas and Fuchssteiner in [31] in connection 
with the Benjamin-Ono Equation. 

Suppose that we have a bi-Hamiltonian system defined by the Poisson tensors 7Ti, tt 2 and 
the Hamiltonians Hi, H 2 . Assume that 7Ti is symplectic. We define the recursion operator 
1Z = 7r 2 7pf 1 , the higher flows 

X, = W- I xi , (47) 

and the higher order Poisson tensors 

7T- = K^TTi . 

For a non-degenerate bi-Hamiltonian system, master symmetries can be generated using a 
method due to W. Oevel [32]. 

Theorem 5 Suppose that X is a conformal symmetry for both tv 2 and Hi, i.e., for some 
scalars X, /i, and v we have 



£x 7ri = A7Ti, Cx K2 = ^2, C>x Hi = vHi 



Then the vector fields 

are master symmetries and we have 
(a) 



x, = n l x. 







(b) 



(c) 



Cx^j = (^ + 0' - * - 2)(/i - A))vr i+i 



[X i ,^] = (^-A)(j-i)X i+j . 



3 A N TODA LATTICE 

3.1 Definition of the system 

Equation (2) is the classical, finite, nonperiodic Toda lattice. This system was investigated 
in [33], [34], [35], [36], [37], [38], [39] and numerous of other papers that are impossible to 
list here. 

This type of Hamiltonian was considered first by Morikazu Toda [39]. The original Toda 
lattice can be viewed as a discrete version of the Korteweg-de Vries equation. It is called 
a lattice as in atomic lattice since interatomic interaction was studied. This system also 
appears in Cosmology. It appears also in the work of Seiberg and Witten on supersymmetric 
Yang-Mills theories and it has applications in analog computing and numerical computation 
of eigenvalues. But the Toda lattice is mainly a theoretical mathematical model which is 
important due to the rich mathematical structure encoded in it. 
Hamilton's equations become 

Qj = Pi 

pj = e qj ~ 1 ~ qj — e qj ~ qj+1 

The system is integrable. One can find a set of independent functions {Hi, . . . , H^} which 
are constants of motion for Hamilton's equations. To determine the constants of motion, 
one uses Flaschka's transformation: 

a % = i e 3(*-*+i) , k = - X - Vi . (48) 

Then 



hi = a { (b i+ i - h) 
h = 2(a 2 i -a 2 i _ 1 ) . 



(49) 



These equations can be written as a Lax pair L — [B,L], where L is the Jacobi matrix 



and 
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This is an example of an isospectral deformation; the entries of L vary over time but the 
eigenvalues remain constant. It follows that the functions Hi = 4-trZ/ are constants of 
motion. We note that 

N 1 

H i = S bi = ~o(Pi + Vi + • • • + Pn) , 
i=i 1 

corresponds to the total momentum and 



N 



N-l 



H 2 = H(q u ...,q N ,Pu---,PN) = -J2 b2 i + J2 a i 



i=l 



is the Hamiltonian. 

Consider H 2N with coordinates (qi, . . . , qN,Pi, ■ ■ ■ ,Pn), the standard symplectic bracket 



N ( 1 

{f,g} s = ( ; 



N ' df_dg_ _ df_dg_ 

dqidpi dpidqi 

and the mapping F : R 2JV — > R 2 ^" 1 defined by 

F : (q 1 ,...,q N ,p 1 ,...,p N ) -> (ai, . . . , ajv-i, h, . . . , b N ) 
There exists a bracket on R 2 ^" 1 which satisfies 

{f,g}oF = {foF,goF} s . 
It is a bracket which (up to a constant multiple) is given by 

{a i} bi} = 
{oi, = ai ; 



(50) 



(51) 



all other brackets are zero. H\ = b\ + b 2 + . . . + b^ is the only Casimir. The Hamiltonian in 
this bracket is H 2 — \ tr L 2 . We also have involution of invariants, {Hi,Hj} = 0. The Lie 
algebraic interpretation of this bracket can be found in [2]. We denote this bracket by tti. 
The quadratic Toda bracket appears in conjunction with isospectral deformations of Ja- 
cobi matrices. First, let A be an eigenvalue of L with normalized eigenvector v. Standard 
perturbation theory shows that 

,2 \T 



VA = (2v 1 v 2 , • . . , 2v N ^v N , v 1 ,...,v N ) := U 

A d\ 
* dai ' ' ' ' ' db?j ■ 

n2 u x = Xm U y 



where VA denotes (4^-, ■ ■ ■ , -ir-)- Some manipulations show that U x satisfies 





— !<2jaj + i 


{ai, k} 


= —Oik 


{a t , b i+ i} 


= aib i+1 


{bi,b i+1 } 


= 2a?; 



where tt\ and 7r 2 are skew-symmetric matrices. It turns out that tx\ is the matrix of coefficients 
of the Poisson tensor (51), and 7r 2 , whose coefficients are quadratic functions of the a's and 
6's, can be used to define a new Poisson tensor. The quadratic Toda bracket appeared in 
a paper of Adler [40] in 1979. It is a Poisson bracket in which the Hamiltonian vector field 
generated by H\ is the same as the Hamiltonian vector field generated by H 2 with respect 
to the 7Ti bracket. The defining relations are 



(52) 



all other brackets are zero. This bracket has det L as Casimir and Hi — tr L is the Hamil- 
tonian. The eigenvalues of L are still in involution. Furthermore, tt 2 is compatible with tti. 
We also have 

7T 2 ViZ, = 7riViZ, + i . (53) 

These relations are similar to the Lenard relations for the KdV equation; they are generally 
called the Lenard relations. Taking / = 1 in (53), we conclude that the Toda lattice is 
bi-Hamiltonian. In fact, using results from [15], we can prove that the Toda lattice is 
multi-Hamiltonian: 

7T 2 V#i = 7Ti VH 2 = 7T V# 3 = 7T-1 Vif 4 = • • • (54) 

Finally, we remark that further manipulations with the Lenard relations for the infinite Toda 
lattice, followed by setting all but finitely many a^, bi equal to zero, yield another Poisson 
bracket, n 3 , which is cubic in the coordinates; see [41]. The defining relations for 7r 3 are 

{ai,ai + i} = aiCii + ibi + i 

{(ii,bi} = -dib\-a\ 

{ai,b i+ i = aib 2 +1 + af , s 

{ai,b i+2 } = aid 2 +l 

{cti + i,bi} = —a 2 ai + i 

{bi,b i+1 } = 2a 2 (bi + b i+1 ) ; 



all other brackets are zero. The bracket tt 3 is compatible with both tt\ and tt 2 and the 
eigenvalues of L are still in involution. The Casimir for this bracket is trL -1 . 



The multi-Hamiltonian structure of the Toda lattice is well-known. The results are usually 
presented either in the natural (q,p) coordinates or in the more convenient Flaschka coor- 
dinates (a, b). In the former case the hierarchy of higher invariants are generated by the 
use of a recursion operator [42], [43]. In the later case one uses master symmetries as in 
[9], [10]. We have to point out that chronologically every result obtained so far was done 
first in Flaschka coordinates (a, b) and then transferred through the inverse of Flaschka's 
transformation to the original (q,p) coordinates. This is to be expected since it is always 
easier to work with sums of polynomials than with sums of exponentials. 
The sequence of Poisson tensors can be extended to form an infinite hierarchy. In order 
to produce the hierarchy of Poisson tensors one uses master symmetries. The first three 
Poisson brackets are precisely the linear, quadratic and cubic brackets we mentioned above. 
If a system is bi-Hamiltonian and one of the brackets is symplectic, one can find a recursion 
operator by inverting the symplectic tensor. The recursion operator is then applied to the 
initial symplectic bracket to produce an infinite sequence. However, in the case of Toda 
lattice (in Flaschka variables (a, b)) both operators are non-invertible and therefore this 
method fails. The absence of a recursion operator for the finite Toda lattice is also mentioned 
in Morosi and Tondo [44] where a Nijenhuis tensor for the infinite Toda lattice is calculated. 
Recursion operators were introduced by Olver [45]. 

3.2 Multi— Hamiltonian structure 

In the case of Toda equations, the master symmetries map invariant functions to other 
invariant functions. Hamiltonian vector fields are also preserved. New Poisson brackets are 
generated by using Lie derivatives in the direction of these vector fields and they satisfy 
interesting deformation relations. We give a summary of the results: 

• There exists a sequence of invariants 

H±, H 2 , H 3 , . . . , 

where — |Tr V , 

• a corresponding sequence of Hamiltonian vector fields 

Xl J X.2 5 X3 5 • • • 5 

where Xi = Xm, 

• a hierarchy of Poisson tensors 

TTl, 7I "2, 7T 3 , . . . , 

where 7Tj is polynomial, homogeneous, of degree i. 

• Finally, one can determine a sequence of master symmetries 

Xi, x 2 , x 3 , . . . , 

which are used to create the hierarchies through Lie derivatives. 
We quote the results from refs. [9], [10]. 



Theorem 6 

i) Kj, j > 1 are all Poisson. 

ii) The functions Hi, i > 1 are in involution with respect to all of the iij. 
hi) XiiHj) = {i+j)H l+J , % > —1, j > 1. 

iv) Lx.-Kj = (j - i- 2)ir i+j , i > -I, j > 1. 

v) [Xi, X j ] = (j-i)X i+j ,i>0,j>0. 

vi) nj V Hi = 7Tj_i V where 7Tj denotes the Poisson matrix of the tensor Tij. 
To define the vector fields X n one considers expressions of the form 

L = [B, L] + L n+1 . (56) 

This equation is similar to a Lax equation, but in this case the eigenvalues satisfy A = \ n+1 
instead of A = 0. 

We give an outline of the construction of the vector fields X n . We define X_i to be 



N f\ 

v^vTr^x;-. 

and X to be the Euler vector field 

N-l Q N q 

We want X± to satisfy 

Xi(Tr L n ) = nTr L n+1 . 
One way to find such a vector field is by considering the equation 

L = [B, L] + L 2 . 



(57) 



(58) 



(59) 



(60) 



Note that the left hand side of this equation is a tridiagonal matrix while the right hand 
side is pentadiagonal. We look for B as a tridiagonal matrix 



B 



(li Pi \ 

"i 72 fa 

a 2 73 fa ■■ 

w ; ••. ••• •••/ 



(61) 



We want to choose the a i: fa and % so that the right hand side of equation (60) becomes 
tridiagonal. One simple solution is a n = — (n + l)a n , fa — (n + l)a n , 7„ = 0. The vector 
field X 1 is defined by the right hand side of (60): 



JV-l 







N 







Xl= ^ dn da- + ^ bn ^r> 



n=l 



n=l 



dbr, 



(62) 



where 



a n = -na n b n + (n + 2)a n b n+1 (63) 

b n = (2n + 3)a* + (1 - 2n)a£_ 1 + 6* . ( 64 ) 
To construct the vector field X 2 we consider the equation 

L = [B, L] + L 3 . 

The calculations are similar to those for X\. The matrix B is now pentadiagonal and the 
system of equations slightly more complicated. The result is a vector field 

N-l Q N _ Q 

^ da n 56 n 

where 

d n = (2-n)a^„ 1 a„+(l-n)a„6^+a n b n 6„ + i + (n+l)a„a^ +1 + (n+l)a n 6^ +1 +a^+o-„a n (6 n+ i-6 n ) 

b n = 2(7„a^-2 ( 7 n _ 1 a2_ 1 +(2n+2)a^ n +(2n+l)a^ n+1 ++(3-2n)a2_ 1 fe n _ 1 +(4-2n)a2_ 1 fe n +^ , 
with 

n-l 
i=l 

and o"i = 0. 

We continue the sequence of master symmetries for n > 3 by 

[X 1 ,X n _ 1 ] = (n-2)X n . (65) 

3.3 Properties of X n and 7r n . 

It is well known that ix\, 7r 2 , 7r 3 satisfy Lenard relations 

7r n Vif, = 7r n _iVff,+i , n = 2, 3 V/ . (66) 

We want to show that these relations hold for all values of n. We denote the Hamiltonian 
vector field of Hi with respect to the nth bracket by In other words, 

X? = [*»,#!]• (67) 
We prove the Lenard relations in an equivalent form. 



Proposition 2 = 



Proof: To prove this we need the identity 

[x 1 , x r] = (n-3) X r l + (/ + i)xr + i (68) 

which follows easily from X\{Hj) — (I + l)Hi +i and (67). Therefore, 
(n-3) X r +1 = [^i,xT]-(i + l)xnn 

= [Xi,xr+i 1 -a+i)xr + i 

= (n - 4) X r +1 + (/ + 2)xS - (I + l)xr+i (69) 
= (n - 4) X r +1 + (/ + 2) X r + i - (I + l)xT+i 

= (^-3)xF + i 

Using the Lenard relations we can show that the functions H n are in involution with respect 
to all of the brackets 7r n . 

Proposition 3 {Hi,Hj} n = 0, where { , } n is the bracket corresponding to ir n . 
Proof: First we consider the Lie-Poisson Toda bracket. We have 

{H 1 ,H j } = Vj, (70) 
since Hi is a Casimir for Suppose that Hj} = V j. 

t{H t ,H 3 } = {X^Hi-^Hj} 

= -Ix^XuHi^]] 

= [X 1 ,{^ i _ 1 ,fr j -}] + [^_ 1 ,C7 + l)xJ + i] (71) 

= (j + l){^_ 1; ff j+1 } 

= . 

Now we use induction on n. Suppose 

{H h H j } n = Q (72) 

{ifj,i?j} n+ i = xr +1 (-^j) 

= Xi+l(Hj) 

— {Hi + l,Hj} n 

= . 

Of course, one can prove the involution of integrals without using master symmetries. We 
present the classical proof: 

In this proof the symbol ir n stands for the bundle map 7r* defined by (18). We first prove the 
involution of constants in the Lie-Poisson Toda bracket. We use the basic Lenard relation 

7T2 dH\ = 7Ti dH.2 ■ 

Since Hi is a Casimir in the linear bracket, we have n\dH\ = 0. We calculate 

{H i ,H j } 1 = < dHi, -KxdHj > 

= — < dHj, TTidHi > 

= — < dHj, -K 2 dHi^i > 

= < df/j-i, 7i2dHj > 

= < dH^i, TV\dHj + i > 

= {Hi, H j+i _ 1 } 1 = . 



It is also easy to show involution in the second quadratic bracket: 

{Hi,Hj} 2 = < dHi, n 2 dHj > 
= < dHi, iridHj +1 > 
= {Hi, H j+1 } 1 = . 

The general result follows from the Lenard relations 7T n dHj = iv n -idHj + i and induction: 

{Hi,Hj} n = <dH i ,n n dH j > 

= < dEi,n n ^ x dE j+1 > 
= {Hi, H j+1 } n _i = . 

■ 

It is straightforward to verify that the mapping 

f(ai, ajv-i, h,...,b N ) = (a x , a N - U 1 + h, . . . , 1 + b N ) (74) 

is a Poisson map between 7r 2 and %i + 7r 2 . Since / is a diffeomorphism, we have the isomor- 
phism 

7T 2 = 7Ti + 7T 2 . (75) 

In other words, the tensor 7r 2 encodes sufficient information for both the linear and quadratic 
Toda brackets. An easy induction generalizes this result: i.e., 



Proposition 4 



n — 1 



J 



*» = E ■ ( 76 ) 



The function tr L? n , which is well-defined on the open set det L ^ 0, is a Casimir for n n , 
for n > 3. The proof uses the Lenard type relation 

7T„VA = AtT^VA (77) 

satisfied by the eigenvalues of L. To prove the last equation, one uses the relation 

TTn E ^ VA * = E A l VA * • ( 78 ) 

But 

X! Air^TrfcVAfc - A fc 7r n _iVA fc ) = , (79) 

for I = 1, 2, . . . , N + 1, has only the trivial solution because the (Vandermonde) coefficient 
determinant is nonzero. 

Proposition 5 For n > 2, trL 2 ~ n is a Casimir for % n on the open dense set detL ^ 0. 
Proof: For n — 3, 

TraVtrL- 1 = 7r 3 £ fc -^VA fc 

/7T„V7 \ , 

(80) 



i 

Efc -^2-A fe 7r 2 VA fc 



= -Efc^iVAfc 
= - 7 r 1 VtrL = xi = • 



For n > 3 the induction step is as follows: 



= E fe (2 - n)l\ n k - VnVAfc 

= Efc(2-n)3^TA fc 7r n _iVA fc (81) 

= ^Elvr^iVtrL 3 -™ 

= . 



3.4 The Faybusovich-Gekhtman approach 

In [46] Faybusovich and Gekhtman find another method of generating the multi-Hamiltonian 
structure for the Toda lattice. Their method is important and will certainly have applications 
to other integrable systems, both finite and infinite dimensional, solvable by the inverse 
spectral transform. Their work shows that the Hamiltonian formalism is built into the 
spectral theory. 

In the case of Toda lattice, the key ingredient is the Moser map which takes the (a, b) 
phase space of tridiagonal Jacobi matrices to a new space of variables (Aj,rj) where Aj is 
an eigenvalue of the Jacobi matrix and rj is the residue of rational functions that appear 
in the solution of Toda equations. The Poisson brackets of Theorem 6 project onto some 
rational brackets in the space of Weyl functions and in particular, the Lie-Poisson bracket 
7Ti corresponds to the Atiyah-Hitchin bracket [47]. We briefly describe the construction: 
Moser in [48] introduced the resolvent 

R(\) = (\I -L)- 1 , 

and defined the Weyl function 

/(A) = (R(X) ei , ei) , 

where e\ = (1, 0, . . . , 0). 

The function /(A) has a simple pole at A = Aj and positive residue at A, equal to rf 

/(A) = £ 



i A — Aj 

The variables (a, b) may be expressed as rational functions of A, and r« using a continued 
fraction expansion of /(A) which dates back to Stieltjes. Since the computation of the 
continued fraction from the partial fraction expansion is a rational process the solution is 
expressed as a rational function of the variables (Aj, r^). The idea of Faybusovich and 
Gehtman is to construct a sequence of Poisson brackets on the space (Aj,rj) whose image 
under the inverse spectral transform are the brackets 7Tj defined in Theorem 6. The Lie- 
Poisson bracket n 1 corresponds to the Atiyah-Hitchin bracket on Weyl functions which in 
coordinate free form is written as 

A rational function of the form is determined uniquely by the dinstinct eigenvalues of 

p(A), Ai, . . . , X n and values of q at these roots. The residue is equal to j^pj and therefore 
we may choose 



X 1 ,...,X n ,q(X 1 ),...,q(X n ) 



as global coordinates on the space of rational functions (of the form - with p having simple 
roots and q,p coprime). We have to remark that the image of the Moser map is a much 
larger set. 

The kth Poisson bracket is defined by 

{A,, ? (Ai)} = -A??(A0 
{q(X t ), q(X 3 )} = {X t , X j } = 0. 

Let us denote this bracket by w^. 

On the other hand in [9], page 108, there is a definition of vector fields on the space of 
eigenvalues of the Jacobi matrix which are projections of the master symmetries X^. They 
are defined by 

N Q 

One verifies easily that these vector fields satisfy the usual Virasoro type relation 

[ej, ej\ = (j — i)e i+ j . 

If we denote by F the function which sends the Jacobi matrix to its eigenvalues then 

dF(X 1 ) = ei 
dF(X 2 ) = e 2 . 

Therefore, it follows by induction that 

dF{X % ) = e t . 

Faybusovich and Gekhtman used the brackets w,i and the vector fields ej to obtain the 
analogue of Theorem 6 in the space of rational functions. The relations obtained correspond 
to the relations of Theorem 6 under the inverse of the Moser map. 

The explicit formulas for for the brackets Wk can be deduced easily from the formulas in [46]. 
They are 

\fii \}k = XiTi 

{A^ Xj}k = . 

In a recent paper [49] Vaninsky has also explicit formulas in (Aj, r*) coordinates for the initial 
bracket w±. 



3.5 A Theorem of Petalidou 

Finally, we mention an interesting result of Petalidou [50]. She proves the following Theorem: 
Suppose that (M, Ao,Ai) is a bi-Hamiltonian manifold of odd dimension and let p be a 
point in M of corank 1. If there exists locally an infinitesimal automorphism Zq of Aq which 



is transverse to the symplectic leaf through p and a vector field Zi which depends on a 
parameter t such that 

[A 1 ,Z 1 ] + [Z 1 ,^]AZ 1 = , 

and 

[A ,Z 1 ] + [A 1 ,Z ] + [Z 1 ,-]AZ = , 

then one can find a symplectic realization of both A and A 1 by a pair of symplectic brackets 
A , Ai given by 

d 

Ai = A i + Z l A— i = l, 2 . 
ot 

This result applies in the case of the Toda lattice by taking Zq = X_\ and Z\ = Xq. 
She obtains symplectic realizations of 7Ti and ty 2 . Furthermore, she constructs symplectic 
realizations of the whole sequence 

7Tl,7T2, 7T3, • • ■ 

of Theorem 6. The corresponding symplectic sequence is given by 

d 

The tensors 7r fc may be generated by a recursion operator since the initial tensor, which is a 
multiple of A , is invertible. 

3.6 A recursive process of Kosmann Schwarzbach and Magri 

In [51] Y. Kosmann-Schwarzbach and F. Magri consider the relationship between Lax and 
bi-Hamiltonian formulations of integrable systems. They introduce an equation, called the 
Lax-Nijenhuis equation, relating the Lax matrix with the bi-Hamiltonian pair and they show 
that every operator that satisfies that equation satisfies also the Lenard recursion relations. 
They derive the multi-Hamiltonian structure of the Toda lattice by defining a matrix M 
and a vector Ao which arise by manipulating the Lax-Nijenhuis equation. They show that 

7T 2 = M?Ti + X ® A , 

where X is the Hamiltonian vector field X2- In the next step of the recursive process they 
show that 

7T 3 = MtY 2 + X ® Al , 

where Ai = MA . In general, 

7T l+1 =MTi i +X®M {l ' 1) \ . 

4 LIE GROUP SYMMETRIES OF THE TODA LATTICE 

Sophus Lie introduced his theory of continuous groups in order to study symmetry properties 
of differential equations. His approach allowed a unification of existing methods for solving 
ordinary differential equations as well as classifications of symmetry groups of partial and 
ordinary differential equations. A symmetry group of a system of differential equations is a 



Lie group acting on the space of independent and dependent variables in such a way that 
solutions are mapped into other solutions. Knowing the symmetry group allows one to 
determine some special types of solutions invariant under a subgroup of the full symmetry 
group, and in some cases one can solve the equations completely. Lie's methods have been 
developed into powerful tools for examining differential equations through group analysis. 
In many cases, symmetry groups are the only known means for finding concrete solutions 
to complicated equations. The method applies of course to the case of Hamiltonian or 
Lagrangian systems, both autonomous and time dependent. Recently, the immense amount 
of computations needed for determining symmetry groups of concrete systems has been 
greatly reduced by the implementation of computer algebra packages for symmetry analysis 
of differential equations. The symmetry approach to solving differential equations can be 
found, for example, in the books of Olver [52], Bluman and Kumei [53], Ovsiannikov [54] 
and Ibragimov [55]. 

Some properties of master symmetries are clear: They preserve constants of motion, Hamil- 
tonian vector fields and they generate a hierarchy of Poisson brackets. We are interested in 
the following problem : Can one find a symmetry group of the system whose infinitesimal 
generator is a given master symmetry? In the case of Toda equations the answer is negative. 
However, in this section we find a sequence consisting of time dependent evolution vector 
fields whose time independent part is a master symmetry. Each master symmetry X n can 
be written in the form Y n + tZ n where Y n is a time dependent symmetry and Z n is a time 
independent Hamiltonian symmetry (i.e. a Hamiltonian vector field). 

In other words, we find an infinite sequence of evolution vector fields that are symmetries 
of equations (49). We do not know if every symmetry of Toda equations is included in this 
sequence. 

We begin by writing equations (49) in the form 



that generate the symmetry group of the Toda equations. The first prolongation of v is 




2aJ_! = . 



We look for symmetries of Toda equations, i.e. vector fields of 



the form 




d_ 

,j w 3 



N-l q N 

pr ( 1 )v = v+ j — + J2g j 

j=l Ua 3 j=l 




pr«(r,) = pr (1) (A,) = 0. 



Therefore we obtain the equations 



4>j ~ ™jOi+i - bj) + (j>j{bj - bj+i) + atfj - a^+i = , 



(82) 



ipj - 2f{a) - a^_ x ) - Aarfj + 4a 7 _i0 j _i = . (83) 



We first give some obvious solutions : 

i) r = 0, cj)j = 0, ipj = 1. This is the vector field 



m^) r = — 1, <j)j = 0, ipj — 0. The resulting vector field is the time translation — whose 
evolutionary representative is 

N ~\ d " d 

j=i j=i 

This is the Hamiltonian vector field xh 2 - It generates a Hamiltonian symmetry group. 
Hi) t = — 1, <pj = a,j, ipj = by Then 

v = + g a ^ + g V = ~dt + x ° • 

This vector field generates the same symmetry as the evolutionary vector field 

Xo + txn 2 ■ 

We next look for some non obvious solutions. The vector field X\ is not a symmetry, so we 
add a term which depends on time. We try 

4>j = -jajbj + (j + 2)ajb j+1 + t( aj a 2 j+1 + a$ +1 - a)_ x a,j - aft]) 

ipj = (2j + 3)a 2 + (1 - 2j)aJ_ 1 + b) + t(2a$b j+1 + 2a) - 2aJ_ 1 a j - 2a 2 j _ 1 b j ) , 
and r = 0. 

A tedious but straightforward calculation shows that ipj, ipj satisfy (82) and (83). It is also 
straightforward to check that the vector field 

^ 0J Da--^Ob 3 

is precisely equal to X 1 + txn 3 ■ The pattern suggests that X n + txu n+2 is a symmetry of 
To da equations. 

Theorem 7 The vector fields X n + tx n +2 ar & symmetries of Toda equations for n > — 1. 

Proof : Note that xh x = because Hi is a Casimir for the Lie-Poisson bracket. We use the 
formula 

[X n , X i] = (l-l)Xn+i ■ (84) 
In particular, for / = 2, we have [X n , % 2 ] = Xn+2- 



Since the Toda flow is Hamiltonian, generated by X2, to show that Y n = X n + tx n +2 are 
symmetries of Toda equations we must verify the equation 



dt 

But 



BY 

" n + [ X 2,Y n } = 0. (85) 



lfr + \X2,Y n ] = + [ X 2, X n + tXn+2] 



Xn+2 - [X n , X2] (86) 

= Xn+2 — Xn+2 = . 



5 THE TODA LATTICE IN NATURAL COORDINATES 



In this section we define the positive and negative Toda hierarchies for the Toda lattice in 
(q,p) variables. We follow reference [15]. 



5.1 The Das-Okubo-Fernandes approach 

Another approach, which explains the relations of Theorem 6 is adopted in Das and Okubo 
[42] , and Fernandes [43] . In principle, their method is general and may work for other finite 
dimensional systems as well. This approach was also used in [56] by da Costa and Marie 
in the case of the Relativistic Toda lattice. The procedure is the following: One defines 
a second Poisson bracket in the space of canonical variables (q 1 , . . . ,q N ,pi, . . . ,p N ). This 
gives rise to a recursion operator. The presence of a conformal symmetry as defined in Oevel 
[32] allows one, by using the recursion operator, to generate an infinite sequence of master 
symmetries. These, in turn, project to the space of the new variables (a, b) to produce a 
sequence of master symmetries in the reduced space. 
Let Ji be the symplectic bracket (50) with Poisson matrix 



Ji = 



I 
-I 



where I is the N x N identity matrix. We use J\ = 4J X . With this convention the bracket J\ 
is mapped precisely onto the bracket 7Ti under the Flaschka transformation (48). We define 
J 2 to be the tensor 

A B 
-B C 



J: 



2 



where A is the skew-symmetric matrix defined by for % < j, B is the diagonal 

matrix (—pi, — p 2 , • • • , — Pn) and C is the skew-symmetric matrix whose non-zero terms are 
Ci,i+i — —Ci+i,i — e qi ~ qi+1 for i — 1, 2, . . . , N — 1. We define J 2 = 2J 2 . With this convention 
the bracket J 2 is mapped precisely onto the bracket 7r 2 under the Flaschka transformation. 
It is easy to see that we have a bi-Hamiltonian pair. We define 

h x = -2(pi +p 2 + ■ ■ ■ +p N ) , 

and hi to be the Hamiltonian: 

N i JV-1 

i=i z i=\ 

Under Flaschka's transformation (48), h\ is mapped onto 4(6i + 6 2 + . . . + b^) = 4 trL = 4H\ 
and hi is mapped onto 2 trL 2 = 4iJ 2 . Using the relationship 

7T 2 V Hi = TTiV H 2 , 

which follows from Proposition 2, we obtain, after multiplication by 4, the following pair: 

JiV h 2 = J 2 V hi . 
We define the recursion operator as follows: 

11 = J2J1 1 ■ 



The matrix form of 7Z is quite simple: 

n= 1 -( B ~ A ) 

2\C B J 

This operator raises degrees and we therefore call it the positive Toda operator 
coordinates, the symbol x% is a shorthand for xk- It is generated as usual by 

X l = w-'xi ■ 

In a similar fashion we obtain the higher order Poisson tensors 

J. = TV- 1 Ji . 
We finally define the conformal symmetry 

N Q N Q 

^o = g(iV-2 i + 1 )^ + g P ,^. 
It is straightforward to verify that 

£-z Ji = —<h , 
£ Z() J2 = . 

In fact, Z is Hamiltonian in the J2 bracket with Hamiltonian function ^Y^iLiQi] see [43]. 
This observation will be generalized in 5.3. 
In addition, 

Z (hi} = hi 
Z (h 2 ) = 2h 2 . 

Consequently, Z is a conformal symmetry for J\, J 2 and h\. The constants appearing in 
Oevel's Theorem are A = — 1, /i = and v — 1. Therefore, we end-up with the following 
deformation relations: 

[Zi,hj] = (i+j)h i+j 

L ZiJj = (j ~ i ~ 2)Ji+j 

\ z h z i\ = U -i)Zi+j ■ 
Switching to Flaschka coordinates, we obtain relations iii) — v) of Theorem 6. 



(87) 
• In (q,p) 



5.2 The negative Toda hierarchy 

To define the negative Toda hierarchy we use the inverse of the positive recursion operator 
1Z. We define 

J\f = 7Zr l = J x J 2 l . 

Obviously we can use the same conformal symmetry Z = K and take A = 0, /i — — 1 and 
v = 2. In other words the role of A and n is reversed. We define the vector fields 

K i =N^K =Af i Z z = 1,2, . . . 



example, Y_i — Ki — J\fZ = —2Y^^ =1 This vector field, in (a, b) coordinates, is given 



which are master symmetries. We use the convention Y_i = K, for i = 0,1,2,.... For 

a 

hv 

N Q 

X_ 1 =V// 1 = VtrL = ^ — 



This is precisely the same vector field (57). In that section, X-i was constructed through a 
different method. Similarly, the vector field Z corresponds to the Euler vector field (58): 

N-l Q N Q 

da, ^ dh 

Note: We use the symbol Y, for a vector field in (q,p) coordinates and X, for the same vector 
field in (a, b) coordinates. Similarly, we denote by J, a Poisson tensor in (p, q) coordinates 
and 7Tj the corresponding Poisson tensor in (a, b) coordinates. The index i ranges over all 
integers. 

We now calculate, using Oevel's Theorem: 

[Y_ h Y_ 3 ] = [K it Kj] = (jm - A)(j - i)K i+j = (-l)(j - i)K i+J = (i - j)Y_ {i+j) . 
Letting m = —i and n = —j we obtain the relationship 

[Y m ,Y n ] = {n-m)Y m+n , (88) 
for all m, n negative. The same relation holds in Flaschka coordinates. In other words 

[X m , X n ] = (n - m)X m+n Vm,neZ- . 

This last relation may be modified to hold for any two arbitrary integers m, n. We suppose, 
without loss of generality, that j > i and consider the bracket of two master symmetries 
Ki = Y_i and Zj = Yj, one in the negative hierarchy and the second in the positive hierarchy, 
i.e. 

Kl = N l Z = R- l Z , 

and 

Zj = n j z . 

We proceed as in the proof of Oevel's Theorem (see [43]): First we note that 
C Zo n = (C Zo J 2 ) Jf 1 - J 2 J^Czo Ji Jf 1 = (fi- \)K ■ 



On the other hand 



c Zi M = ^ Zo (JiJ 2 1 ) = (\-^W ■ 



Finally, 



[Y-i, Y j \ = [K i ,Z j ] = \Af i Z ,WZ ] 
= M'Czq (TZ j ) Z - K j £ Zo (Af l ) Z 

AFjifi - \)n j z - nn{\ - ll)N 1 z 

= - A)7#"% - i(X - /x)7#-% 



= {i + 3)( i i-\)W- i Z Q 



In the case of Toda lattice \i = and A = — 1, therefore 



\Y- i ,Y j \ = (i + j)Y j 



We deduce that (88) holds for any integer value of the index. 

We define Wi = J 3 _j. This is necessary since the conclusions of Oevel's Theorem assume 
that the index begins at % — 1 and is positive. We compute 



£y_, J-j = C Kl W J+3 = {n+(j+3-2-i){n-\))W i+j+3 = (i-j-2)W i+j+3 = (i-j-2)J_ {i+j) . 



for n, m negative integers. Switching to Flaschka coordinates we deduce that the relation 
iv) of Theorem 6 holds also for negative values of the index. In other words 



Again, a straightforward modification of the proof of Oevel's Theorem shows that the last 
relationship holds for any integer value of m, n. We have shown that conclusions iv) and v) 
of Theorem 6 hold for integer values of the index. In fact, it is not difficult to demonstrate 
all the other parts of Theorem 6. 

Theorem 8 The conclusions of Theorem 6 hold for any integer value of the index. 
Proof. 

We need to prove parts i), ii), Hi) and vi) of the Theorem. 

i) The fact that J n are Poisson for n E Z follows from properties of the recursion operator. 
The similar result in (a, b) coordinates follows easily from properties of the Schouten bracket, 
and the fact that J n and n n are F— related. We have n n = F*J n , therefore 



The vanishing of the Schouten bracket is equivalent to the Poisson property. 

Hi) The case where i and j are both of the same sign was already proved. We next note 

that X n (X) = X n+1 if A is an eigenvalue of L. This follows from equation (56) which is used 



Letting m = —i and n = —j we obtain 




L x^j = {j-i- 2)ir i+j , i < 0, j < . 



K, 7r n ] = [F*(J n ),F m (J n )] = F*[J n , J n ] = F,(0) = . 



to define the vector fields X n for n > 0. We would like to extend the formula X n (X) = A n+1 
for n < 0. Since X^i(X) = 1 we consider X_ 2 . We look at the equation 

[X_ 2 ,X n ] = (n + 2)X n _ 2 . 

We act on A with both sides of the equation and let X_ 2 (A) = /(A). We obtain the equation 

(n + l)A/(A)-/'(A)A 2 = (n + 2) . 

This is a linear first order ordinary differential equation with general solution 

/(A) = I + cA" +1 . 

Since n is arbitrary, we obtain f(X) — j. In order to calculate X_ 3 (A) we use 

X-3 = — [X-i, X-2] ■ 

We obtain 

X_ 3 (A) = X_ 2 X^{\) 

The result follows by induction. 
Finally we calculate 

Uh 3 ) = -.x, (£ aj) = y(E ^ A i) = E ^"^(Afc) = E Aj- 1 Aj+ 1 = E K +j = (i+m 

vi) First we note that ttj V Hi = Kj-i Vffj+i, holds for i, j of the same sign. More 
generally, in the positive (or the negative) hierarchy we have the Lenard relations (77) for 
the eigenvalues, i.e. 

tvj V A, = AiTTj-iV A, . (89) 
Assume now that % < 0, j > 0. The calculation is straightforward: 

^V±£A| = EArSVA, 
= EA^-iVAfc 
= 7r i _ 1 V i £ I EAS+ 1 . 

Therefore, 

TTjVJJi = TTj-iV^i+i • (90) 

In the case z > and j < we use exactly the same calculation but use (89) for the negative 
hierarchy. 

ii^ It is clearly enough to show the involution of the eigenvalues of L since Hi are functions 
of the eigenvalues. We prove involution of eigenvalues by using the Lenard relations (90). 
We give the proof for the case of the bracket ttj with j > but if j < the proof is identical. 
First we show that the eigenvalues are in involution with respect to the bracket Let A 
and \i be two distinct eigenvalues and let U, V be the gradients of A and /i respectively. We 
use the notation { , } to denote the bracket 7Ti and ( , ) the standard inner product. The 
Lenard relations (89) translate into 7r 2 U = \niU and n 2 V = /jtti V. Therefore, 



-X_xX_ 2 (A)= -*_!(-) = -. 



{A,//} = (mU,V) = j{n 2 U,V) 

= ~(U,ir 2 V) = ~{U, f jLir 1 V) 
= -^(U,7r 1 V) = ^(n 1 U,V) 

= 

Therefore, {A,/x} = 0. To show the involution with respect to all brackets Hj, and in view 
of part iv) of Theorem 6, it is enough to show the following: Let fa, fa be two functions in 
involution with respect to the Poisson bracket n, let X be a vector field such that X(fa) = ff 
for % — 1,2. Define a Poisson bracket w bj w = Cxtt- Then the functions fa, fa remain in 
involution with respect to the bracket w. The proof follows trivially if we write w = Cx^ in 
Poisson form: 

{fa, fa} w = X{fa, fa} n - {fa, X(fa)} n - {X(fa), fa} n . 

■ 

Remark: We should point out that 

H n = -trL n , 

n 



makes sense for n ^ but it is undefined for n — 0. The reader should interpret the 
formulas involving H as a degenerate case, i.e. H = ^f- = ^ = oo. Therefore the result 
of X- n (H n ) = N where N is the size of L. It makes sense to define 



X m (if ) = lim-X m (trL") . 

n-»0 n 



For example, X_i(H ) is calculated by X_i ^ tr L n ^j = tr L n 1 . Taking the limit as n — > 
gives X_i{H Q ) = trL -1 = —H_i which is the correct answer. 



5.3 Master integrals and master symmetries 

In this section we prove some further results and give some specific examples. 
In subsection 5.1 we noticed that Zq is Hamiltonian with respect to the J 2 bracket with 
Hamiltonian function f — \ Hf=i Qi- This observation is due to Fernandes [43]. This type of 
function is called a master integral. It is not a constant of motion, but its derivative is. We 
generalize the result as follows: 

Theorem 9 The master symmetry Y n , n G Z is the Hamiltonian vector field of f with 
respect to the J n+2 bracket. 

Proof. We will prove the result for the positive hierarchy Z n = Y n but the proof for Y_ n = K n 
is similar. As a first step we show that 

Z n (f) = V n > . 



We recall that 



N f) N f) 



Since 



N 



£(iV+l-2i)=0, 



i=l 



we obtain 

i N 1 AT 

W) = ^o(Eft) = 9(E^o(ft)) = o. 
z i=i z i=i 

By examining the form (87) of the recursion operator 1Z we deduce easily that the % com- 
ponent of Z 1 is 



^i(ft) = ~2 



In other words, the vector 
is the product where 



(^-2i + l)pi + Eft"-Eft- 
(Zi(gi), . . .,Zi(gjv)) 



/^o(?i) 1 1 
-1 1 

-1 -1 Z (g 3 ) 

V -1 



1 \ 



1 

-1 Z (q N )J 



and P is the column vector (pi,P2, • • • ,PnY- Note that X^=i a y — and X^=i % — ~Z (q, 
Therefore, 



1 1 

z i(f) = ^( z i{qi) + ■■■ + Z^qif)) = - E a«Pj ,, 

In the same fashion one proves that ^(/) = 0. 
For n > 2, we proceed by induction. 



i 



n-2 



[Z U Z., 



n-l 



Therefore, 



Z n (f) 



' [Zx, Z n _x] (/) = — ^ (ZxZ^x/ - Z^iZx/) = , 



n-2 



n-2 



by the induction hypothesis. 

To complete the proof of the Theorem, it is enough to show 



Z n — [Jn+2, f] , 



where [ , ] denotes the Schouten bracket. 
First we note that 



[[J n+1 , /], Z x \ + [[/, Z x \, J n+1 ] + [[Z u J n+1 ],f] = 
due to the super Jacobi identity for the Schouten bracket. Since 

[Z 1 ,f] = Z 1 (f) = 0, 
the middle term in the last identity is zero. We obtain 

[Z 1 ,[J r n. 1 ,f]] = [[Z 1 ,J n+1 ],f] . 

Finally, we calculate using induction: 

Z n = Z n -l] 

71/ Zi 

-[Z^Jn+lJ]] 



n-T 
= ^[[Z 1) ^ 1 ],/] 

' ([(n-2)J n+2 ,/]) 



n-2 

= [Jn+2, f] ■ 

■ 

The result of the Theorem is striking. It shows that the master symmetries are determined 
once the Poisson hierarchy is constructed. Of course one requires knowledge of the function 
/. The function / may be constructed by using Noether's Theorem: The symmetries of the 
Toda lattice in (q,p) coordinates have been constructed in [57], at least for two degrees of 
freedom. The Lie algebra for the potential of the Toda lattice with TV degrees of freedom is 
five dimensional with generators 



X 1 



d_ 

at 

N Q 

i=2 u i % 

N Q 
N Q 



The non-zero bracket relations satisfied by the generators are 

[X U X 2 ] = X 1 

[X\ , X 5 ] = X4 

[X2, X 3 ] = —2X 4 

[X 2 ,X 5 ] = x 5 

1X31X4] — —2X 4 

IX3, X 5 ] = —2X 5 . 



This Lie algebra L is solvable with L« = [L, L] = {X u X 4 , X 5 }, = {X 4 } and = {0}. 
We examine the symmetry X 5 . 

A corresponding time dependent integral produced from Noether's Theorem is 

1 N i N i 
1 = o^ qi ~ 2 f ^ Pi = ^ + 1 tkl ' 

i=l i=l 

Motivated by the results of [58], [59], it makes sense to consider the time independent part of 
/ which is precisely the function /. It is an interesting question whether this procedure works 
for other integrable systems as well. We also remark that the integrals are also determined 
from the knowledge of the Poisson brackets and the function /. For example, it follows easily 
from Theorem 9 that ^ 

hi+i = . {K, f}3 , 
where { , } 3 denotes the cubic To da bracket. 



5.4 Noether Symmetries 

We recall that Noether's Theorem states that for a first order Lagrangian, the action inte- 
gral J t * 2 Ldt is invariant under the infinitesimal transformation generated by the differential 
operator, known as a Noether symmetry, 

d N d 

* = T£l + EQi£- (92) 
dt fr[ dqi 

if there exists a function F, known as a gauge term, such that 

fir N Fit N Fit 

F = Tj^ + Qijr- + E (4 - ftT) P + TL. (93) 
dt S dqi fr[ V / dqi 

When the corresponding Euler-Lagrange equation is taken into account, equation (93) can 
be manipulated to yield the first integral 



F — 



N 



PiT 

TL + E (Qi - QiT) j- 



(94) 



Thus to every Noether symmetry there is an associated first integral. Consider the La- 
grangian to be of the form 

1 N 

L = -E^ 2 -^(<7i>92,---,gjv) . (95) 
1 i=i 

We summarize the results of [60] in the following Theorem: 

Theorem 10 Let X be the N x 1 vector with entries Qi, x the vector with entries qi and 
b the vector with entries bi(t). Let A be an N x N skew- symmetric matrix with constant 
entries. Finally we denote by In the N x N identity matrix. If X , given by (92), is a Noether 
symmetry then the infinitesimals must be of the form 

T = T(t) 

x = ( A + if lN ) x+b < 96 > 



and the gauge term is restricted to 



„ ld 2 T^ 2 "dbAt) 
The associated first integral I is equal to 

N 

F + TH - QiPi , 
i=i 

where H is the Hamiltonian. 

By examining the form (91) of the generators for the Toda lattice we conclude, using Theorem 
10, that only Xi, X4 and X5 are Noether symmetries. The corresponding integrals provided 
by Noether's Theorem are the Hamiltonian H = h 2 , the total momentum hi — p± + . . . + Pn 
and / + th 1 . 

In order to obtain more integrals we consider generalized Noether symmetries. That is, 
the infinitesimals in (92) do not just depend on t, qi, . . . , qw but also on qi, . . . , q^. For 
Lagrangians of the form (95) for one, two and three degrees of freedom, all the possible 
point Noether symmetries are classified in [60]. The following results are from [61]. 
In the case of generalized symmetries, we can take without loss of generality T — 0. Hence, 
we consider operators of the form 

G = |>£ < 98 > 

where the infinitesimals of (92) and (98) are related by 

Vi = Qi- Qi T ■ 
Using (95) and (98), Noether's condition (93) becomes 

N N N N ( N N \ 

ft + Yl + 2 9<4 = 2 ViKi + ft [vit + Yl QjViqj + S ■ (") 

1=1 1=1 1=1 1=1 y j=i j=i j 

We consider equation (99) in the case of the Toda lattice with two degrees of freedom. 
By assuming various forms of the r]i (i.e. linear, quadratic or arbitrary) we can solve this 
equation and produce the following integrals one of which (J 3 ) is new: 



92 



h=Pi+P2 , h = (pi - P2) 2 + 4e qi - 



P1-P2 + \fh ( rrqi + q2 



Pi-P2-Vh V P1+P2J 

Note that H = + I 2 ) and that the function G = which appears in the exponent 

of ^3 satisfies {G, H} = 1. 

The existence of the integral J 3 shows that the two degrees of freedom Toda lattice is super- 
integrable with three integrals of motion J 2 , 13}. A Hamiltonian system with iV degrees 
of freedom is called super-integrable if it possesses 2N — 1 independent integrals of motion. 
Of course these integrals cannot be all in involution. Based on this computation for 2 degrees 
of freedom we make the following conjecture: 
Conjecture: The Toda lattice is super-integrable. 



5.5 Rational Poisson brackets 



The rational brackets in (q,p) coordinates are given by complicated expressions that are 
quite hard to write in explicit form. When projected in the space of (a, b) variables they give 
rational brackets whose numerator is polynomial and the denominator is the determinant of 
the Jacobi matrix. We give examples of these brackets and master symmetries for N = 3. 
For example, the tensor J is a homogeneous rational bracket of degree 0. It is defined by 

In the case of three particles the corresponding bracket ir is given as follows: First define 
the skew-symmetric matrix A by 



a 12 = -^aia 2 (^3 + h - b 2 ) 

a 13 = a 1 (a 2 2 -b 2 h) 
an = -a-i(a\ - b-ibz) 

2 

Cil5 — (l\d2 
2 

&23 — ~ a \ a >2 

a 24 = a 2 (al-b 1 b 3 ) 
a 25 = -a^aj-hh) 
a 34 = -2a\b 3 

a 3 5 = 

a 45 = — 2a\b\ . 

The matrix of the tensor 7r is defined by 

= ddi- 4 < 100) 

where detL = bib 2 b 3 — a\bi — a\b 3 . This formula defines a Poisson bracket with one single 
Casimir H 2 = |trL 2 . The bracket is defined on the open dense set detL ^ 0. Taking 
H 3 = |tr L 3 as the Hamiltonian we have another bi-Hamiltonian formulation of the system: 

7Ti dH 2 = 7T dH 3 . 

In fact we have infinite pairs of such formulations since 

tx 2 dH\ = 7Ti dH 2 = 7Tq dH 3 = 7r_i dH/± = . . . 

The explicit formulas for the vector fields Xi and X 2 are given in 3.1 therefore we give an 
example for the vector field X_ 2 . In the case iV = 3 it is given by 



i * a ^ 



where 



n = 501(61-62-263) 

r 2 = 2 a2 (^3 - 26i - 6 2 ) 

si = 6 2 6 3 - a\ — a| 

«2 = 6163 + af + 

S3 = 61 6 2 - a\ - a\ . 

Finally, we consider the Casimirs of these new Poisson brackets. 

Theorem 11 The Casimir of 7i n in the open dense set det L 7^ is Tr L 2 ~ n for all n 7^ 2. 
The Casimir of 7r 2 is det L. 

Proof: 

For n > 1 the result was proved in Proposition 5. Therefore, we only have to show that the 
Casimir of 7T_ m is tr L m+2 (m > 0). This follows from (90) and the fact that ifi = Tr L is 
the Casimir for the Lie- Poisson bracket tt\\ 

= TTiVFl = 7r V# 2 = 7T_iV# 3 = • • • • 



6 GENERALIZED TODA SYSTEMS ASSOCIATED WITH SIM- 
PLE LIE GROUPS 



In this section we consider mechanical systems which generalize the finite, nonperiodic Toda 
lattice. These systems correspond to Dynkin diagrams. They are special cases of (1) where 
the spectrum corresponds to a system of simple roots for a simple Lie algebra. It is well 
known that irreducible root systems classify simple Lie groups. So, in this generalization for 
each simple Lie algebra there exists a mechanical system of Toda type. 
The generalization is obtained from the following simple observation: In terms of the natural 
basis qi of weights, the simple roots of A n _i are 

Ql - Q2, 52 - 93, • • • , Qn-1 ~ Qn ■ (101) 

On the other hand, the potential for the Toda lattice is of the form 

e qi ~ q2 + e q2 ~ q3 + ... + e qn - 1 ~ qn . (102) 

We note that the angle between q^i — and — q i+ i is ^ and the lengths of — q i+1 are 
all equal. The Toda lattice corresponds to a Dynkin diagram of type A n _ x - 

More generally, we consider potentials of the form 

U = Cl e Mq) + ... + cie fl(q) (103) 

where ci, . . . , q are constants, fi{q) is linear and I is the rank of the simple Lie algebra. For 
each Dynkin diagram we construct a Hamiltonian system of Toda type. These systems are 



interesting not only because they are integrable, but also for their fundamental importance 
in the theory of semisimple Lie groups. For example Kostant in [2] shows that the integration 
of these systems and the theory of the finite dimensional representations of semisimple Lie 
groups are equivalent. 

For reference, we give a complete list of the Hamiltonians for each simple Lie algebra. 



C n 



Go 



Ep. 



E 7 



Ea 



1 n 

H = -J2p 2 j+ e qi ~ 92 + ■■■ + e 9 "" 1 - 9 " 

2 i 



1 n 

H = -Y jP 2 j + e qi - q2 + ■■■ + e qn - 1 - q " + e q " 

2 J 



H = - Y,P 2 j + e qi ~ q2 + ■■■ + e 9 - 1 " 9 " + e 2qn 



1 n 

ff — _ p 2 _|_ g9i~92 _|_..._|_ e qn-i~Qn _|_ e q„-i+q n 

2 1 3 



1 3 

H=-J2p 2 + e 9l ~ Q2 + e- 2qi+q2+q3 

2 l 



1 4 

JJ — _ \^ p 2 -4- e 9l _ 92 _|_ g92~93 _|_ g<?3 _|_ g 5(94 -Ql- 92-93) 

2 1 5 



1 8 4 

H — _ + V" g*-9j + l _|_ g-(91+92) _|_ gi(-<?i +92 + . -+97-98) 

2 1 1 



1 8 5 

/J — _ + ^ e 9j— _|_ g"(9l+92) _|_ g|(-9l +92 + . ..+97-98) 

2 1 1 



1 8 6 

H — _ ^2p 2 + V" e 9j— _|_ g-(91+92) _|_ gi(-9l +92 + - -+97-98) 

2 1 1 



We should note that the Hamiltonians in the list are not unique. For example, the A 2 
Hamiltonian is 



H = \v\ + \v\ + \v\ + e 91 " 92 + e" ■ (104) 



An equivalent system is 



H(Q h = -Pi + -Pi + e vf(^i+<W + e" 2 ^ . (105) 
The second Hamiltonian is obtained from the first by using the transformation 



Qi = 


V2, 

— {qi + Q2 -2q 3 ) 


(106) 


<?2 = 


Ve, , 
— \Q2 -qi) 


(107) 


Pi = 




(108) 


P2 = 




(109) 



Another example is the following two systems, both corresponding to a Lie algebra of type 
D 4 : 

4 n 2 

^2 — + e gi + e 92 + e qz + e^ qA ~ qi ~ q2 ~ q ^ 
i=i 2 

4 2 

El. _|_ e 9l-92 _|_ e Q2-Q3 _|_ e 93-<?4 _|_ e 93+94 

i=l 2 

Finally, let us recall the definition of exponents for a semi-simple group G. An excellent refer- 
ence is the book by Collingwood and McGovern [62]. Let G be a connected, complex, simple 
Lie Group G. We form the de Rham cohomology groups H l (G,C) and the corresponding 
Poincare polynomial of G: 

where di =dim H l (G, C). A Theorem of Hopf shows that the cohomology algebra is a finite 
product of I spheres of odd dimension where / is the rank of G.. This Theorem implies that 

p G (t)=ri(i+^ +i ) . 

i=i 

The positive integers {ei, ei) are called the exponents of G. One can also extract 

the exponents from the root space decomposition of G. The connection with the invariant 
polynomials is the following: Let Hi, H 2 , Hi be independent, homogeneous, invariant 
polynomials of degrees mi,m 2 ,..,mj. Then = m, — 1. The exponents of a simple Lie 
group are given in the following list: 

4.-1 



1,2,3, ...,n-l 



1,3,5, ....,2n- 1 



1, 3, 5, 2n — 3, n — 1 



1,5 



1,5,7,11 



1,4,5,7,8,11 



1,5,7,9,11,13,17 



1,7,11,13,17,19,23,29 



7 5 n TODA SYSTEMS 

7.1 A rational bracket for a central extension of 5 n -Toda 

In this subsection we show that the B n Toda system is bi-Hamiltonian by considering a 
central extension of the the corresponding Lie algebra in analogy with gl(n, C) which is a 
central extension of sl(n, C) in the case of A n Toda. 

Another way to describe these generalized Toda systems, is to give a Lax pair representation 
in each case. It can be shown that the equation L = [B, L] is equivalent to the equations of 
motion generated by the Hamiltonian H2 = |tr L 2 on the orbit through L of the coadjoint 
action of £?_ (lower triangular group) on the dual of its Lie algebra, B*_. The space B*_ 
can be identified with the set of symmetric matrices. This situation, which corresponds to 
sl(n, C) = A n _i can be generalized to other semisimple Lie algebras. We use notation and 
definitions from Humphreys [63]. 

Let Q be a semisimple Lie algebra, $ a root system for Q, A = {cci, . . . , a{\ the simple roots, 
h a Cartan subalgebra and Q a the root space of a. We denote by x a a generator of Q a . 
Define 

B- =/i©£& • 

a<0 



G2 

Er 

E 8 



There is an automorphism a of Q, of order 2, satisfying a(x a ) = x_ a and a(x_ a ) = x a . 
Let /C = {x G £/ | o"(cc) = — x}. Then we have a direct sum decomposition Q = £>_ © /C. 
The Toda flow is a coadjoint flow on B*_ and the coadjoint invariant functions on Q*, when 
restricted to B*_ are still in involution. The Jacobi elements are of the form 



L = b i h i + J2 a i + X -*i) ■ 



i=l 



We define 

i 

B ^ ] (-^aj aj) • 

i=l 

The generalized Toda flow takes the Lax pair form: 

L=[B,L] . 

The B n Toda systems were shown to be Bi-Hamiltonian. The second bracket can be found in 
[10]. It turned out to be a rational bracket and it was obtained by using Dirac's constrained 
bracket formula (29). The idea is to use the inclusion of B n into A 2n and to restrict the 
hierarchy of brackets from A 2n to B n via Dirac's bracket. Straightforward restriction does 
not work. We briefly describe the procedure in the case of B 2 . 



The Jacobi matrices for A4 and B 2 are given by 



La 4 



and 



Lb 2 — 



a 1 



V 













\ 


a\ 


b 2 


a 2 











a 2 


h 


a 3 











a 3 




04 


V 










bj 



CLi 

b 2 
a 2 







a 2 

h 

-a 2 







-02 

26 3 - h 
— a\ 







-Oi 

26 3 - h J 



(110) 



(in) 



Note that L^ 4 lies in gl(4, C) instead of sl(4, C). Therefore we have added an additional 
variable in L B2 . We define 



Pi — 0,1 + 04 

p 2 = a 2 + a 3 

P3 = b 1 + b 5 - 2b 3 

Pa = b 2 + b 4 - 2b 3 . 

It is clear that we obtain B 2 from At by setting pi = for i = 1,2,3,4. We calculate the 
matrix P = {pi,Pj}. The bracket used is the quadratic Toda (52) on A4. 



{Pl,P2> 


= ^(aia 2 -a 3 a 4 ) 


{Pi,Ps} 


= (X465 — a\b\ 


{Pi,Pa} 


= a x b 2 — 0464 


{P2,Ps} 


= 2(a 3 b 3 - 2a 2 b 3 ) 


{P2,Pa} 


= a 3 6 4 + 2a 3 b 3 - 2a 2 b 3 - a 2 b 2 


{P3,Pa} 


= -2a\ - Aa\ + Aa\ + 2a\ . 



If we evaluate at a point in B 2 we get 



{Pl,Ps} 

{P2,Ps} 
{P2,P4} 
{^3,P4} 



Therefore the matrix P is given by 



and P 1 is the matrix 



P- 1 



( 


2a\b 3 4a 2 b 3 

\-2a x b 3 Qa 2 b 3 



( 







-2a x b 3 
2a x b 3 
-4a 2 b 3 
-Qa 2 b 3 
. 



-2a x b 3 
-4a 2 6 3 







2a x b 3 \ 
-Qa 2 b 3 


J 



V 



10ai&a 
1 

10a2 63 






3 


1 


\ 


10ai6 3 


5ai&3 





1 


1 




10a2&3 


10a 2 63 




1 










5ai&3 




1 

10a2 &3 








/ 



Using Dirac's formula (29) we obtain a homogeneous quadratic bracket on B 2 given by 

aia 2 (3b 3 — b 2 — 26i) 



{a l} a 2 } 

{ai, h} 

{ai,b 2 } 

{01,63} 
{a 2 ,&i} 



IO63 

-a l (m 1 b 3 -2b 1 b 2 -?>b\- a\) 
IO63" 

a!(106 2 6 3 - 3b 2 2 - 2b x b 2 - Aa 2 2 - aj) 



106, 



01(62 - 61) 



a 2 (2M 3 - 26ib 2 + a\) 
IO63 



{a 2 ,b 2 } 
{a 2 ,6 3 } 
{61,62} 
{61,63} 



-02(86263 - 3bj - 6al - Aa\) 

10^ 

a2(6 3 - 6 2 ) 
5 

10a% - 3a% - 1a\b x - 3a?&i 
56s" 

2a? 
5 



{62,63} = -{a\ - a\) . 

The bracket satisfies the following properties which are analogous to the quadratic A n Toda 
(52). 

i) It is a homogeneous quadratic Poisson bracket. 

ii) It is compatible with the B 2 Lie- Poisson bracket. 

Hi) The functions H n = ^ tr L n are in involution in this bracket. 

iv) We have Lenard type relations 7r 2 Vifi = HiVH i+ i where 7Ti, 7t 2 are the Poisson matrices 
of the linear and quadratic B 2 Toda brackets respectively. 

v) The function det L is the Casimir. 



7.2 A recursion operator for Bogoyavlensky-Toda systems of type B n 

In this section, we show that higher polynomial brackets exist also in the case of B n Toda 
systems. We will prove that these systems possess a recursion operator and we will construct 
an infinite sequence of compatible Poisson brackets in which the constants of motion are in 
involution. 

The Hamiltonian for B n is 



We make a Flaschka-type transformation 



1 n 

H = - Y,p] + e gi ~ q2 + ■■■ + e 9 - 1 ^" + e qn . (112) 



a^lelfe-^i), a n = ie^" (113) 

h = ~Pi ■ 



Then 

ai = a>i (6i+i -h) i = l,...,n 
hi = 2(a?-a?_i) i = l,...,n, 

with the convention that ao = 6 n+ i = 0. 

These equations can be written as a Lax pair L = [B, L], where L is the symmetric matrix 



(114) 



( bi a\ 



a„_i 

a n-l b n 

a n 





-a n 



~b n 



(115) 



V -ai -61 y 

and I? is the skew-symmetric part of L (In the decomposition, lower Borel plus skew- 
symmetric). We note that the determinant of L is zero. 

The mapping F : R 2n — > R 2n , (qi,Pi) — > (oi,6i), defined by (113), transforms the standard 
symplectic bracket into another symplectic bracket 7Ti given (up to a constant multiple) by 



(116) 



{ai,6j} = -aj 
{cii,b i+1 } = ^ . 

It is easy to show by induction that 

1 1 22 2 

det 7Ti = a x a 2 . . . a n . 

The invariant polynomials for B n , which we denote by 

H 2 , H 4 , . . . H 2n 

are defined by H 2 i = ^ Tr L 2 \ The degrees of the first n (independent) polynomials are 
2, 4, . . . , 2n and the exponents of the corresponding Lie group are 1, 2, . . . , 2n — 1. 

We look for a bracket 7r 3 which satisfies 



7T 3 V H 2 = 7Ti V #4 • 

Using trial and error, we end up with the following homogeneous cubic bracket 7r 3 . 

i — 1, 2, . . . , n — 1 



;ii7) 









= — aib\ — af 


{a n ,& n } 


= -a n b 2 n -2al 


{fli, ^+2} 


= dia i+1 


{flj, b i+ i\ 


= dib 2 +l + a? 


{ai, bi-i} 


= -a^Oj 


{bi,b i+l } 


= 2a?(6i + 6 i+ i) 



(118) 



We summarize the properties of this new bracket in the following: 

Theorem 12 TTie bracket 7r 3 satisfies 

1. 7r 3 is Poisson 

2. 7r 3 compatible with tt\ . 

3. H 2i are in involution. 



Define 1Z = . Then 1Z is a recursion operator. We obtain a hierarchy 

7Ti, 7T 3 , 7T 5 , . . . 

consisting of compatible Poisson brackets of odd degree in which the constants of motion are 
in involution. 

4. 7V j+2 V H 2i = TTj V H 2i+2 V i, j . 
The proof of this result is in [10]. 

It is interesting to compute the cubic Poisson bracket in (p, q) coordinates. We will see that in 
the expression for the master symmetry the exponents of the corresponding Lie-group appear 
explicitly. We reproduce the formula for the cubic Poisson bracket in (p, g)-coordinates from 
[11]. 



{q%,Qi- 

{Pi, ft- 

{Pn, Qn-2} = {Pn, q n -3} 



-l} = {ft,ft-2> 

-2} = {Pi,qi-3} 



{Qi,Qi} 
{Pi,Qi} 

= {Pn,Ql} 



2pi i = 2,...,n 

2(e qi ~ 1 ~ qi — e qi ~ qi+1 ) 



3, . . . , n — 1 



Ui,Pi} 
{q n ,p n } 

{Qi+uPi} 

{qi,Pi+i} 

{q n -i,Pn} 

{Phpi+l} 



p\ + 2e qi ~ qi + 1 s 
pi + 2e qn 

e Qi-qi+l 

2^e qi+1 ~ qi+2 — e qi ~ qi+1 ) 



l,...,n- 1 



l,...,n-2 



In (p, g)-coordinates, J\ is the (symplectic) canonical Poisson tensor, h 2 is the Hamiltonian , 
J3 is the cubic Poisson tensor for B n and Z is the conformal symmetry for both Ji, J 3 and 
h 2 . 

So, with 

nan Q 

we have 



£-z Ji — —Jii £>z J3 — <hi £>z h 2 — 2h 2 . 

We obtain a hierarchy of Poisson tensors, master symmetries and invariants which are ob- 
tained using Oevel's Theorem. For example, we have 

[Z i , Xj ] = (2j + l)Xi +j (119) 

and the coefficients of the first n independent Hamiltonian vector fields correspond to the 
exponents of a Lie group of type B n . 



7.3 Bi— Hamiltonian formulation of B n systems 



Following the procedure outlined in the introduction we obtain a bi-Hamiltonian formulation 
of the system. In other words, we define 7r_! = A/Vx = niir^ 1 % 1 and we use it to obtain the 
desired formulation. 

We illustrate with the B 2 Toda system. In this case det 7i\ = a\a\ and det 7r 3 = a\a\ A 2 = 
det 7Ti A 2 where 

A = a\ + 2a\a\ + 2a\b\ + b\b\ - 2a\b x b 2 . 



The explicit formula for 7r_x is 



7T_ 



where 



A 





aia 2 6 2 
ai(&2 + cl\ 1 



2a 2 ) 



—aia 2 b 2 


—afa 2 



V-ai(fe 2 + a 2 + 2a 2 ) a 2 (6? + 2a 2 ) 



-ai (61 + a 2 + 2al) ai (6 2 + a? + 2a|) \ 

a 2 a 2 -a 2 (6 2 + 2a^) 

-2a 2 (6! + 6 2 ) 

2a 2 (^i + 62) / 



This bracket is Poisson by construction. We will prove later that it is compatible with 7Ti. 
We note that A = Vdet 7Z and it is also equal to the product of the non-zero eigenvalues of 
L. Using the rational bracket 7r_i we establish the bi-Hamiltonian nature of the system, i.e. 



mVH 2 = 7r_iV#4 • 



8 C n TODA SYSTEMS 



We now consider C n Toda systems. We will prove that these systems also possess a recursion 
operator and we will construct an infinite sequence of compatible Poisson brackets as in the 
B n case. We also show that the systems are bi-Hamiltonian. 

8.1 A recursion operator for Bogoyavlensky— Toda systems of type C n 

The Hamiltonian for C n is 

1 11 

H=-Y / p 2 i + e qi ~ q2 + ■■■ + e^- 1 -^ + e 2qn . (120) 

2 1 

We make a Flaschka-type transformation 

o i = ^e3(*- 9i + 1 ) , a n = ^=e qn (121) 



The equations in (a, b) coordinates are 



hi = ^ (b i+ i - bi) i = 1, . . . , n - 1 

a n = -2a n b n (122) 
bi = 2(a 2 -a 2 i _ 1 ) i = 1, 



with the convention that a^ = 0. 

These equations can be written as a Lax pair L — [B,L], where L is the matrix 



a n -l 






a n _i b n 


a n 




a n 


-b n 


— a n-l 




—a>n-\ 





V 



-ai 



-ai 

-61 y 



(123) 



and 5 is the skew-symmetric part of L. 

In the new variables Oj, 6j, the canonical bracket on R 2n is transformed into a bracket 7Ti 
which is given by 



{cii,bi} 
{ai, b i+1 } 
{a n , b n ^ 



—ai i — 1, 2, . . . , n — 1 
Oj z = 1, 2, . . . , n — 1 



(124) 



-2a 



n 



The invariant polynomials for C n , which we denote by 

H 2 , H 4 , . . . H 2n , 

are defined by H 2 i = ^ Tr L 2 \ 

We look for a bracket 7r 3 which satisfies 

7T 3 V if 2 = tti V #4 . 

The bracket 7r 3 was obtained in [12]: 



(125) 



{aj, tti + i} 

{a n -i, a n } 
{a i} k} 
{a n , b n } 
{a i} b i+2 ] 
{ai, b i+ i} 
{a n ~i, b n } 
{a i} bi-i} 

{a n , b n -i} 
{bi,b i+1 } 



aia i+ ib i+1 i = 1,2, ..,n - 2 



-2a n & 2 _ 2a 3 



1,2, 



n — 1 



= ^n— 1 a n —\b n a n —\a n 
= -a\_ x ai 



(126) 



2a?(6i + 6 i+ i) . 

We summarize the properties of this bracket in the following: 

Theorem 13 The bracket 7r 3 satisfies 

1. 7r 3 is Poisson 

2. 7r 3 is compatible with tti. 

3. H 2i are in involution. 



Define 1Z = . Then 1Z is a recursion operator. We obtain a hierarchy 

7Ti, 7T 3 , 7T 5 , . . . 

consisting of compatible Poisson brackets of odd degree in which the constants of motion are 
in involution. 

4. n j+2 V H 2i = TTj V H 2i+2 V i, j . 

The proofs are precisely the same as in the case of B n . 

Even though it is not necessary to work in (p, ^-coordinates, we reproduce the formulas 
from [11] for completeness. 



{Qi,Qi-i} = {Qi,Qi-i} = ■ ■ ■ = = 2Pi i = 2,...,n 

Kfc-2} = {Pi,qi-s} = ■■■ = {Pi,Qi} = 2(e*- 1 "« - e*-'^ 1 ) i = 3,...,ra-l 

{Pn, qn-2} = {Pn, qn-s} = ■■■ = {p n , qi} = 2e q -^ q " - Ae 2q " 



(127) 



(128) 



{qi,Pi} = p 2 + 2e"+ 1 i = l,...,n-l 

{q n ,Pn} = p 2 n + 2e 2qn 

{Qi+i,Pi} = e qi ~ 9i+1 

{qi,p i+ i} = 2e q '+ 1 ~ q '+ 2 - e qi ~ qi+1 i = l,...,n-2 

{q n -i,Pn} = 4e 2 «" - e qn - x ~ qn 

{PhPi+i} = -eP- qi+1 (pi+Pi+i) . 

For C n , the conformal symmetry is 

and we have the same constants as in the case of B n : 

£z Jo = —Jo, Cz Ji = Ji, Cz Q Ho = 2Hq. 
The relations of Oevel's Theorem are the same of the B n Toda 

[Z i , Xj ] = (2j + l)Xi +j (129) 
[Zi,Zj] =2(j-i)Z i+j (130) 



C Zi J j ^{2{j-i)-l)J i+j . 
Note that (129) gives a method of generating the exponents. 



(131) 



8.2 Bi Hamiltonian formulation of C n systems 



In order to show that the C n Toda systems are bi-Hamiltonian we define 7r_i = ttiTT^tti. 
This is the second bracket required to obtain a bi-Hamiltonian pair. We illustrate with a 
small dimensional example, namely C 3 . The explicit formula for 7r_i is the following: Let A 
be the skew-symmetric 6x6 matrix defined by the following terms: 

a 12 = aia 2 (b 2 al + bjb 2 - a\b 3 - b 3 b\ + a\b 2 + b\b 3 ) 

a 13 = -2aia 3 (a|fo 2 + bjb 3 - 6 3 6|) 

a i4 = ai{b\a\ + a\ + b\h\ — 2a\b 2 b 3 + a\ a\ + a\b\) 

ai5 = — ai(2a 2 — 2a 2 b 2 b 3 + a\a\ + a\b\ + a\b\ + b\b\) 

«i6 = ai(2a 2 + 2a§&^ + a 2 a\ — 2a 2 b 2 b 3 — a 2 b\ — 2b\a 3 ) 

a 23 = 2a 2 (b\ + a\)b 3 a 3 

a 24 = -a\a 2 (a 2 3 + b\ - 2b 2 b 3 + a\ - 2b x b 3 ) 

«25 = a 2 {2a\a\ + 2a\b\ — 2a\b 2 b 3 + 2a\a\ — 2a\bib 3 + a\b\ + b\b\ + a\b\) 

a 26 = -a 2 (2ajal + 2a\b\ + a\b\ + 2a\a\ + a\ - 2a\b x b 2 + b\b\) 

a 34 = 2a\a 3 (a\ — 2bib 3 — 2b 2 b 3 ) 

a 35 = — 2a 3 {2a\a\ — 2a\bib 3 + a\b\ — 2a\b 2 b 3 ) 

a 36 = 2a 3 {2a\b\ + 2a|a? + af - 2a^i6 2 + ^6?) 

a 45 = 2af (a§6i + 6 2 a| + b x b\ + £>§6 2 ) 

a 46 = 2a\(a 2 bi — 2a|6i — 2b 2 a\ — a 2 b 3 ) 

a 56 = 2{2a\a\bi + 2a\b 2 a\ + 2a\a 2 b 3 — 2a\a 2 bi + a 2 6^6 3 + a 2 b\b 2 ) . 
The Poisson tensor 7r_x is of the form 

= delL A ' 

where 

det L = VdetTZ = 2a 2 2 b\ b 2 b 3 -2a\a\ b x b 3 -a 2 3 b\ b 2 2 +2a 2 3 b x b 2 a\-a\al-b\ b\ b\+2a\b x b 2 b 2 3 -a\ b\-a\b\ 



As in the case of B 2 we have 



mVH 2 = 7r_iV#4 • 



9 D n TODA SYSTEMS 

In this section, we show that higher polynomial brackets exist also in the case of D n 
Bogoyavlensky-Toda systems. Using Flaschka coordinates, we will prove that these sys- 
tems possess a recursion operator and we will construct an infinite sequence of compatible 
Poisson brackets in which the constants of motion are in involution. We also show that the 
system is bi-Hamiltonian. 

9.1 A recursion operator for D n Bogoyavlensky-Toda systems in Flaschka co- 
ordinates 

The Hamiltonian for D n is 

1 n 

H= J2p 2 j + e qi ~ q2 + ■ ■ ■ + e^- 1 '^ + e qn - 1+qn n>4. (132) 

2 1 



We make a Flaschka-type transformation, F : R 2n — > R 2ri defined by 

F : (gi, . . . ,q n ,pi, . . . ,p n ) -> (ai, . . . , a„, h, . . . , b n ) , 

with 



0+ = - e ^ qi ~ qi+l) 
2 



i = 1,2,..., n- 1, a n = i e lfe-i+^) 



(133) 



Then 



&n-l 



-a n (fe n _i + 6 n ) 
2 (a 2 - a 2 _ x ) 
2(a 2 + a 2 _! - a 2 _ 2 ) 



z = 1, 2, . . . , n . 

i — 1, 2, . . . , n — 1 
i — 1,2, . . . ,n — 2 and i = n 



(134) 



These equations can be written as a Lax pair L = [B, L], where L is the symmetric matrix 



Oi 





O n -l 


-On 












a n 


-On 







— o n -i 





a n 


— On-l 





V 



—a! 
-Oi -0i / 



(135) 



and i? is the skew-symmetric part of L (In the decomposition, lower Borel plus skew- 
symmetric) . 

The mapping F : R 2n — > R 2n , (qi,Pi) — > (aj,6j), defined by (133), transforms the standard 
symplectic bracket J into another symplectic bracket 7Ti given (up to a constant multiple) 
by ' 

{cLi,bi} = -\di i = l,2,...,n 

{ai,6i + i} = \a,i i := 1,2, . . . ,n - 1 (136) 
{a n ,b n -i} = — 2 a n- 

We obtain a hierarchy of invariant polynomials, which we denote by 

H 2 , H A , . . . , H 2n , . . . 

defined by H 2i = \- Tr L 2t . The degrees of the first n — 1 (independent) polynomials are 
2, 4, . . . , 2n - 2. We also define 

P n = VdetL . 

The degree of P n is n. The set {i^2, ^4, H2n-2,Pn} corresponds to the Cheval- 

ley invariants for a Lie group of type D n . The exponents of the Lie group is the set 



{1,3,5,. ..,2n — 3, n — 1} which is obtained by subtracting 1 from the degrees of the invariant 
polynomials. A conjecture of Flaschka states that the degrees of the Poisson brackets is in 
one-to-one correspondence with the exponents of the corresponding Lie group. 
Taking H 2 = |Tr I? as the Hamiltonian we have that 

gives precisely equations (134). In this section, we find a bracket %_i which satisfies 

TT 1 VH 2 = 7T_iV#4 • 

First, we define a bracket 7r 3 which satisfies 

7T 3 V H 2 = TTi V H A , (137) 

and whose non-zero terms are 

{ai,a i+ i} = diCti+ibi+i i = 1,2, . . . ,n - 2 

{a n - 2 ,a n } = a n - 2 a n b n -i 

{di,bi} = -ai(bl + af) i = 1, 2, . . . , n - 2 

{On-lA-l} = -0-n-l{a 2 n-l + + K-l) 

{a n ,b n } = -a n (a 2 n + b\ - al_ ± ) 

{ai,b i+1 } = ai(a 2 + 6- +1 ) i = l,2, ...n-2 

{a n _i,6 n } = a^^a 2 .! + b\ - a 2 ) (138) 

{ai,b i+2 } = aiCL 2 i+1 i = l,2,...,ra-3 

{a„_2,&n} = a„_ 2 (a 2 t „ 1 - a 2 ) 

{ai,6j_i} = -a^cii i = 2,3, . . . ,n - 1 

{a„,6„_i} = -a n (3a 2 _ x + a 2 + 6 2 _ x ) 

{6i,6 i+ i} = 2a 2 (bi + b i+1 ) i = l,2,...,n-2 

{6„_i, 6„} = 2a 2 _ 1 (6 n _i + 6„) + 2a 2 (6 n - 6 n _i) . 

This bracket appeared recently in [13]. We summarize the properties of this new bracket in 
the following: 

Theorem 14 The bracket 7r 3 satisfies 

1. 7T3 is Poisson. 

2. 7r 3 is compatible with tt\ . 

Define 1Z = 7r 3 7rf 1 . Then 1Z is a recursion operator. We obtain a hierarchy 

^1,^3,^5, ■ ■ ■ 

consisting of compatible Poisson brackets of odd degree in which the constants of motion are 
in involution. 

3. All the H 2i and P n are in involution with respect to all the brackets 7Ti, 7r 3 , 7T 5 , 

4- n j+ 2 V H 2i = iTj V H 2i+2 V i, j . 

The proof of 1. is a straightforward verification of the Jacobi identity. We will see later, in 
the next subsection, that 7r 3 is the Lie derivative of 7Ti in the direction of a master 
symmetry and this fact makes 7Ti, 7r 3 compatible. 4- follows from properties of the recursion 
operator. 3. is a consequence of 4 ( see f° r example Proposition 3 for a method of proof). 
The only part which is not obvious is the involution of P n with H n which will be proved at 
the end of next subsection using master symmetries. 



9.2 Master symmetries 

We would like to make some observations concerning master symmetries. 

Due to the presence of a recursion operator, we will use the approach of Oevel. We define 

Z to be the Euler vector field 

~i dcii dbi 
We define the master symmetries Zj by: 

Zj = 7VZq . 

For obvious reasons we use the notations 

X 2i = Zi ,hi = H 2 i+2 , Ifj = 7r 2 i+i , ipi = X2i+2 , i — 0, 1, 2, . . . 

where Xn denotes the Hamiltonian vector field generated by H 2i , with respect to tv-\_. This 
notation is convenient since X 2 is a master symmetry which raises the degrees of invariants 
and Poisson tensors by 2 each time. One calculates easily that 

C Zo U = -n , Czjli = III , £z ho = 2h . 

Therefore Z is a conformal symmetry for Il , IT! and h . The constants appearing in Oevel's 
Theorem are A = — 1 , fx — 1 and z/ = 2. Therefore we obtain 

[Zi, ^] = (1 + 2j) ^ [X 2i , X2 J+ 2] = (1 + 2j) X2(i+i+i) (139) 

[Z, Z,] =2(j- i) [^2i, * 2 ;] = 20- i) X 2( , +J) (140) 

£ Zi (n,) = (2j - 2i - 1) n i+j ^ C X , M (7r 2j+1 ) = (2j -2t-l) 7r 2(i+j)+1 (141) 

Z (fy) = (2 + 2? + 2j) h i+j X 2l (H 2j ) = 2{i+ j) H 2{i+j) . (142) 

Remark 1: The relation (141) implies that £x 2 ( 7I "i) — — 3vr 3 and therefore 7r 3 is Lie-derivative 
of 7Ti in the direction of a master symmetry. This makes 7Ti compatible with 7r 3 (see Lemma 

!)• 

Remark 2: The relation (139) gives a procedure for generating almost all the exponents. As 
we mentioned in the introduction, the last exponent is generated by the application of the 
conformal symmetry on the Hamiltonian vector field corresponding to the Pfaffian of the 
Jacobi matrix. 

It is interesting to note that one can obtain the master symmetry X 2 by using the matrix 
equation 

L = [B, L] + L 3 , (143) 
where L is the Lax matrix (135) and B is the skew-symmetric matrix defined as follows: 



B 






Xi 


y\ 






















x 2 


2/2 














— yi 


— X2 


n 
u 






n 








n 

u 











x n-2 


?/n-2 


y n -i 














— X n -2 





1 


x n 


















%n— 1 










-yn-i 










-y n -i 










%n—l 


-yn-2 























— ^n-2 














y n -i 


y n -2 


^n-2 






























where 



y2 




x 2 

yi 



-y2 
o 

Xi 



X ,r 



yi 

Vn-l 



fi-1 

a i l'^b j + (i + l-n) (hi + b i+1 ) 

n-2 

(i + l—n) aia,i + i, % — 1, 2, . . . ,n — 2 

a n-2a n - 



1,2,. ..,71-1 





-yi 



It is interesting to note that the is a constant times a product of aja^ where the product 
is determined from the Dynkin diagram of a Lie algebra of type D n . 

The matrix B was chosen in such a way that both sides of (143) have the same form. The 
components of the vector field X 2 are defined by the right hand side of (143). 
Finally we note the action of the first master symmetry on 

X 2 (P n ) = P n H 2 . 

Remark: This last result should be expected since the eigenvalues of L satisfy A = A 3 under 
(143). Therefore, 



X 2 (P n ) = X 2 (v / detL) 



X 2 
2 



Ai . . . A n 

l - (Ai . . . A n )~^ (AiA 2 . . . X n + AiA 2 . . . A n + . . . + Ai 
2y/d tL (^^ 2 • • • ^« + ^iA 3 , . . . A n + . . . + Ai . . . A 3 



■•An) 



detL 
2 v / detl 



(xl + xl + ... + x n ) 



= PnH 2 . 



We conclude this subsection by proving the involution of Hi with P n . It is clearly enough 
to show the involution of the eigenvalues of L since P n and Hi are both functions of the 
eigenvalues. It is well-known that the eigenvalues are in involution with respect to the 
symplectic bracket 7Ti. We will give here a proof based on the Lenard relations (137). Let A 
and \x be two distinct eigenvalues and let U, V be the gradients of A and fj, respectively. We 
use the notation { , } to denote the bracket 7Ti and ( , ) the standard inner product. The 
Lenard relations (137) translate into ir 3 U = A 2 7Ti U and 7r 3 V = /i 2 ttx V. Therefore, 



{A,//} = (n^V) 

= ^(^U,V) 

= -^(u,*iV) 

= ^{A,^}. 

Therefore, {A, /j,} = 0. To show the involution with respect to all brackets n 2j+1 and in view 
of (141) it is enough to show the following: Let f\, f 2 be two functions in involution with 
respect to the Poisson bracket 7r, let X be a vector field such that X(fi) = ff for i — 1,2. 
Define a Poisson bracket to by w = jC-xK- Then the functions f\, f 2 remain in involution 
with respect to the bracket w. The proof follows trivially if we write w = Cx^ in Poisson 
form 

{h, f2} w = x{f 1} f 2 } w - {A, x(f 2 )} n - {x(/0, / 2 } w . 



Remark: We have to point out that unlike the case of B n and C n the cubic bracket (138) 
was discovered not by manipulating the left hand side of (137) but through the use of the 
master symmetry X 2 . In other words, we construct the master symmetry X 2 using (143) 
and then compute 7r 3 = — \ i Cx 2 ' K i- 

9.3 A recursion operator for D n Toda systems in natural (q,p) coordinates 

We now define a bi-Hamiltonian formulation for D n Bogoyavlensky-Toda systems in natural 
(qi,Pi) coordinates. This bracket is simply the pull-back of 7r 3 under the Flaschka trans- 
formation (133). After some tedious calculation, we obtain the following bracket in (qi,Pi) 
coordinates: 



{Qi,Qj} = ~2pj , i<3 

{ft, Pi} = ^ + 2e«-*+ 1 i=l,2,...,ra-2 

{g„_i,p„_i} = #U + 2e«"- 1 -«" + 2e«"- 1+ «» 

{g i} Pi-i} = e*- 1 "* i = 2,3,...,ra-l 

{?n,Pn-i} = e'- 1 "** - e«»- 1+ «» 

= _ e «-9i+i + 2e 9!+1 - 9! + 2 i = l,2,...,w-3 

{gn-2,Pn-i} = -e 9 ™- 2-9 ™- 1 + 2e 9 ™- 1 - ffe + 2e <fe - 1+,? ™ 

{?n-i,Pn} = -e qn - 1 ~ q " + e qn - 1+qn 

{qi,Pj} = -2e«- 1- « + 2e«-*'+ 1 l<i<j-l<ra-3 

{g i} p n _i} = -2c 9 ™- 2 - 9 ™- 1 + 2e 9 ™- 1 - 9 ™ + 2e q "- 1+qn i = 1, 2, . . . , n - 3 

{&,Pn} = ^e 9 "- 1- ** + 2e 9 ™- 1+9 " i = 1, 2, . . . , n - 2 

{Pi,Pi+i} = -eF- qi+1 (pi+Pi+i) i = l,2,...,ra-2 

{p n -l,Pn} = -(Pn-1 +Pn)e' 7n - 1 -' 7n + -p„)e 9 ™- 1+< ?" ; 



and all other brackets are zero. 

Denote this Poisson tensor by J\ and let Jo be the standard symplectic bracket. A simple 
computation leads to the following: 

Theorem 15 The bracket J\ satisfies 

1. J i is Poisson. 

2. J\ is compatible with J . 

3. The mapping F given by (133) is a Poisson mapping between J\ and the cubic bracket 7r 3 . 

Thus, in (q,p) coordinates we also have a non-degenerate pair (J , J±) for D n Bogoyavlensky- 
Toda and therefore we may define a recursion operator 1Z = JiJq 1 . We obtain a hierarchy 
of mutually compatible Poisson tensors defined by Jj = TZ 1 Jq. 
The vector field 

z ° = t*k + P {n - i >k- (145) 

is a conformal symmetry for the Poisson tensors Jq and J\ and for the Hamiltonian 

1 n 

h = -Y jP 2 + e qi - q2 + ■■■ + e qn - 1 - qn + e qn - 1+qn . (146) 

2 i 

We compute 

£z Jo = -Jo , £z a h = Ji , £z o h = 2h . (147) 

So Oevel's Theorem applies. With Z { = 1Z l Z , xo — Xh and Xi = ^Xo one calculates easily 
that 

(a) [Zi,Xj] = ( y l + 2j) X i+j 

(b) [Zi,Zj] = 2(j-i)Z i+j 

(c) C Zi (Jj) = (2j-2i-l)J i+J . 

Note that (a) gives the exponents (except one) for a Lie group of type D n . 



The action of the first master symmetry on P n is the same as in Flaschka coordinates: 

Z 1 (P n ) = h P n . (148) 



Finally, we calculate that 
producing the last exponent. 



[Z ,XP n } = (n-l) X P n , (149) 



9.4 Bi-Hamiltonian formulation of Bogoyavlensky— Toda systems of type D n 

In order to show that the D n Toda systems are bi-Hamiltonian we use the same procedure 
as in the previous two cases. The tensors 7Ti and 7r 3 are both invertible and we define 
7r_! = 1 7r 1 . This is the second bracket required to obtain a bi-Hamiltonian pair. We 
illustrate with a small dimensional example, namely -D 4 . The explicit formula for 7r_4 is the 
following: 

Let A be a skew-symmetric 8x8 matrix given by the following terms: 

«i2 = -aia 2 {a\a\b A + 626364 + 2b 4 b 3 b 2 al - 2b 4 b 3 b 2 al + b 2 a\ + a\b 2 - 2b 2 a\a\ + 0^6264 
—636364 + 646^3 — 646^4 + 646364 + 640464 — 64O364 — 036364 — 0^0464) 

fl i3 = _ fl i a 3(626364 — 646^3 — a\b 2 b\ + 646^04 — 646364 — b\a 2 A b 4 + b\a\b 4 — a 2 b 4 + 2a\b 2 b 3 b 4 
—afalbi + a\b 2 a\ — a\a\b 2 — b\b\b 4 + 036363 — b 3 b\a\ + b\a\b 4 + 646364 — 64O363 + 6^6304) 



Oi 4 = 0404(036264 — 636364 + 646303 — 646304 + 646364 + 64O464 — 6^364 — 0^64 + 2O2626364 — O4O264 

+026204 — a\a\b 2 — b\b\b± + 036363 — b 3 b\a\ + b\a\b 4 + b 2 b 3 b4 — 6^0363 + 6^6304) 



fl i5 = —01(204636364 + 264O2O362 + 6^4 + 0^64 — 2b 3 b 2 a\b\ + 6^3 — 2b 3 b 4 b\a\ — 2b\a\a\ — 2b i a\a\b 2 

+b\b 2 4 b\ + 046364 — 204636403 + 204636404 + 0^04 — 2a\a\a\ + a^Og) 



016 = 04(204636404 + 264636404 — 204636403 — 2640^0462 + 2640^362 + 20364 — 26362O264 

+ 0404 + 0403 + 6^4 + 64 O3 - 20^0403 - 26^0^4 + 04636^ - 264O36364 + 646364) 



Oi 7 = 04(26304 - 264636404 + 26^3 + 204636364 - 20364 - 462O4O3 + 26362O264 - 263646303 

-26^ - 26^ + 46? o^ + 2b\a%b A + 6?o 2 ! 6;j - a? o^) 



018 = a\{a\a\a 2 4 — o^OgOg — 26^636404 + 26304 — 2640^0462 — 262O3 + 204636364 — 2640^0362 

2„2 o^2„4 , o^2„4 ql2„2i, r, . r2„2„2 J,2„2„2\ 



+263646^03 - 26? 04 + 26?Og - 26^036364 + 6?020g - 6^0304) 



«23 = 0203(610364 - 616364 - 6^464 - 630^64 - 0^0464 + 0^364 + 6^6364 - 6^363 + 6^6304 

+6364^ + 0463^ — 63CJ1O3 — 646162 + 26461O162 — a\a\b A — 64O1) 

024 = —0204(616364 + 6^464 — 61O364 + 63O164 + 0^0464 — O1O364 + 616364 — 6^0363 + 616304 

+636401 + 046301 — 630^03 — 646163 + 264610162 — a\a\b^ — 64O1) 

025 = —0102(2646303 — 6364 — 26 4 o4 2 6 3 — 04 + 20403 — 03 + 2616364 — 26461O3 + 264O461 — 0364 

+263626^ + 26 4 6 2 o^ - 26 4 o^6 2 ) 

026 — —02(401636404 + 2616364O4 — 2016462O4 — 20^636264 + 201640362 — 261010464 + 261010364 

-26163O164 - 401636403 + 2fli04 + 2aiOg + 61O4 + 61 03 - 4O1O4O3 - 261O3O4 + 2016364 - 

2b\a%b 4 + 6?o^6^ + 2o?a 2 ! 6 4 ! + 6?6 3 : 6 4 ! ) 

027 = 02(0^64 + 2016364O4 + 261636404 — 2016364O3 + 20^04 + 20^03 + 26^4 + 26 2 03 — Aa\a 2 4 a\ 

-461O3O4 - 261036364 + 610364 + 2010364 - 261O16264 + 616364) 

028 = —02(046163 — 2016364O4 — 2616364O4 — 2016364O3 — 20^04 + 20^03 — 26i O4 + 261O3 
— 261O36364 — 2bia\a\b 2 — b\a\b\ + b\a\a\ — b\a\a\ — 0^03 + 0^4 + 261010362) 

O34 = —2030464(6103 + 6162 — 26iOi6 2 + O1O2 + Oi) 



035 = —0301(2646103 — 2616364 — 2640461 + 0^64 — 2636264 — 26462O4 + 2640362 + 2616403 + 
2646364 — 261O363 + 2046163 — 2036364 + O2O3 — 0^04 + 2626364 + 2046263 — 2036263) 

036 = —03(261640362 — 461O1O264 + 2a\a\a\ — 2a\a\a\ + 2016463O4 + 201636364 — 201O46263 
+201036263 + 203616364 + 201610363 + 401036364 — 201636364 — 201046163 — 201616364 — 20^640363 

+261O1O464 - 2610^0364 + 26163O16 2 - 6^026 2 - 201036^ - 61O2O3 + 6^204) 

037 = 03(261640363 — 2bia\a\b^ + 2a\a\a\ — 0^64 — 2a\a\a\ + 203616364 + 201036364 + a\b\b\ 
-2b\a\b\ - 2a\a\b\ + 2b l a\b 2 b\ - 2b x a\a\b 2 - b\a\b\ - 26? o^ + 26?a 2 ! o 4 ! - 6?6^6 2 - 0^3 

+O1O4 + 261010363) 



^38 = 03(2610^263 + 2a\a\a\ + 2a\a\a\ + 0^63 + b\b\b\ + 3a\b\b\ — Qbia^a^ + b^a^bl + b\a\ 
+2blalal + 2b\a\a\ — 2b\a\bzb 2 — 2bia\b 2 b\ + a\a\ + 30^04 — 2bia\a\b 2 ) 

045 = 040^(2616364 — 26461 03 + 264O461 — O264 + 2636264 + 26462O4 — 2640362 + 26164O2 + 2616364 — 

2640363 + 2046163 — 20^6364 + OgOg — 0^04 + 2626364 + 2046263 — 2036263) 

46 = a A {2b\b A a\b 2 — Abia\a\b A + 2(^0^^ — 2a\a\a\ — 2a\b A b 2 a\ — 2a\b^b 2 b\ — 2a\a 2 4 b 2 b ?) + 2a\a\b 2 b^ 

+2026^6364 + 20^610363 + 40^036364 — 20^626364 — 20^046163 — 20^616364 + 20^640362 
-2610^0464 + 2b x a\a\b± - 2bib 3 a(bl + b\a\b\ + 2a{a\h\ - 6?aija§ + 6^0204) 

047 = —04(26^640362 — 2610^264 + 2a\a\a\ + 0^64 — 2a\a\a\ + 20^6^6364 + 20^0^6364 + a\b\b\ 

+2b\a 2 2 bl + 201O2&4 - 2b 1 a\b 2 b 2 A - 2b x a 2 ^a\b 2 - b\a\b\ - 26^203 + 26^0304 + 6^264 

— a^Og + 0^04 + 2610^0362) 

048 = —04(2610^0363 + 2a\a\a\ + 2a\a\a\ + 0^63 + b\b\b\ + a\b\b\ — 2b\a\a\b 2 + 36^0363 

+6^2 + 2b\a\a\ + 2b\a\a\ — 2b\at 2 b- i b 2 — 2b\a?[b 2 b\ + 30^03 + 0^04 — 6610^0362) 

56 = -20^(616364 - 261640363 + 261640463 + 61 04 - 261O3O4 + 61 03 + 626364 + 264636204 - 264636203 

+O462 - 2620403 + 6203) 

057 = —20^(4620403 — 26104 — 26103 — 20462 — O26364 + 2646362O3 — 2646362O4 — 262O3 — 

261640463 + 4610304 + 26164O363 + 610364) 

58 = — 2o^(26i03 — 26104 — 20462 — 264636203 — 264636204 + 26 2 Og + 0^0364 — 26164 O4 63 — 

261640363 — 61O2O4 + 0^0464 + 610^03) 

67 = —Aa\b 2 a\ — 4o^Og6 2 — Ab^a\a\b\ + 80^036204 — Aa\b\a\ — 4a^6iOg + Aa\b\a\b\ — 2a\b\b^b\ 

—2a\b\b\b 2 + 80^610304 + 40^61640363 — 40^61640463 — 463640^6204 + 463640^0362 

068 = —463640^0362 — 463640^6204 — 40^61640363 — 40^61640463 — Aa\b 2 a\ + Aa\a\b 2 — Aa\b\a\ 

+4a^6iOg + Aa\a\a\b A + Aa^a^a^b^ + 4b\a\a\a\ — 4t>ia\a\a 2 4 + 2b\a\a 2 4 b^ 

+26^20364 + 2b\a\b 2 a\ - 26^036203 



a-is = 2a 2 A b^a\ — 2a 1 a 2 b4 — Aa\a\a\b& — Aa\a\o? A b± + 46101036263 — Ab\a\a\b2bz + Ab\bia\a\b2 
—4bia\a\a\ + Abia\a\a\ — 2b\a\b%b\ — 26^646^03 — 2b\b A b\a\ + Abib A a\a\b2 — Ab\a\a 2 A b A 
+2bjb 3 blal - Ab\a\a\b± - Abja^aj + Ab\a\b 2 a\ - 2036301 - 2a\a\b A . 

The Poisson tensor of 7r_i is given by the formula 

det L — P 2 — (036164 + 036162 — 046162 — 61626364 — a 2 a 2 + a 2 a 2 + a 2 6 3 6 4 ) 2 . 
We note that 

det tt 3 = det 7Ti(det L) 2 = P 4 4 , 

and therefore 

det ft = (detL) 2 . 

This formula (as well as the formulas in the previous two cases) indicates that the eigenvalues 
of the recursion operator should be the squares of the non-zero eigenvalues of the Jacobi 
matrix. This is known to hold in the case of the classical A n Toda lattice, see [26]. For a 
general recursion operator, the relation between its eigenvalues and the eigenvalues of the 
Lax matrix has not been fully investigated yet. 
As in the case of B 2 and C 3 we have 

7TiV# 2 = 7T_iVf/ 4 • 

10 Conclusion 

10.1 Summary of results 

The classical, finite, non-periodic Toda lattice is known to be bi-Hamiltonian. Moreover, 
(53) is a multi-Hamiltonian formulation of the system. We have indicated how to obtain 
similar results for the other classical Lie algebras and we have illustrated with some small 
dimensional examples. These examples may be generalized: 

Theorem 16 The B n , C n and D n Toda systems are bi-Hamiltonian. In fact, they are 
multi-Hamiltonian. In each case we define 

M = TTi^ 1 , 

where ~K\ is the Lie-Poisson bracket, 7r 3 is the cubic Poisson bracket and 

7T-(2i-i) =A/Vi i = 1,2,... . 
Then all the brackets are mutually compatible, Poisson and satisfy 



7TiV#2 = 7T_iV#4 = 7T_ 3 V# 6 = 



(150) 



Proof: The proof of (150) is trivial. They are just the Lenard relations for the negative 
hierarchy. The brackets 7r_i, 7t_ 3 , 7r_ 5 , . . . are all Poisson since they are generated by the neg- 
ative recursion operator, Af, applied to the initial Poisson bracket t\\. To prove compatibility 
of all Poisson brackets appearing in (150) we take two brackets ix t and 7r s where t, s are odd 
integers, with t < s < 1. Using condition (b) of Oevel's theorem (for the negative operator) 
we can express irt as the Lie derivative of n s in the direction of a master symmetry. It is 
therefore enough to prove the following simple general result: If ir and a are both Poisson 
tensors and a = Cx^ for some vector field X, then 7r and a are compatible. The one line 
proof uses the super- Jacobi identity for the Scouten bracket: 

[n,a] = [n,[n,X]] = -^[X, [7r,7r]] = 0. 

■ 

We remark that Oevel's theorem applies in all three cases (and for both hierarchies) since 
the Euler vector field X Q (58) is a conformal symmetry for 7Ti, 7t 3 and H 2 . Furthermore, the 
compatibility condition holds for all brackets in both hierarchies. If 7T t and tt s are both in the 
positive hierarchy then the argument of the theorem still works using the positive recursion 
operator. The remaining case, when one of the tensors has negative index and the other 
one positive, can also be proved in a straightforward manner using similar arguments, i.e., 
properties of the Schouten bracket and the fact that the formulas in Oevel's theorem hold 
for any integer value of the index (Theorem 8). The tensor in the positive hierarchy is the 
Lie derivative of the tensor in the negative hierarchy in the direction of a suitable master 
symmetry and the argument of the theorem shows that they are compatible. Therefore, we 
have a more general result: Any two brackets in either the positive or negative hierarchy are 
compatible. 

Remark The compatibility condition follows also from a general result of bi-Hamiltonian 
geometry: If tx and a are two compatible Poisson tensors and n is invertible, then Af = o-k^ 1 
is a recursion operator (i.e. its torsion vanishes) and all the tensors Af l 7T, with i > 0, are 
Poisson and compatible. Using this result one can prove compatibility of all brackets in 
both hierarchies. For example, to show that 7r 5 is compatible with 7r_ 3 we use the fact that 
1Z = 7r 3 7rf 1 = 7r_! (/T_ 3 ) _1 to obtain 7r 5 = 7£ 4 7r_ 3 . Since n_i and 7r_ 3 are compatible and 
Poisson, then 1Z generates a chain of compatible Poisson tensors. In particular n 5 and 7r_ 3 
are compatible. 

10.2 Open problems 

We conclude with some open problems and some possible directions of research for systems 
related to the Toda lattice. 

• Exceptional Toda lattices 

The case of exceptional simple Lie groups is still an open problem. It is also a much 
more difficult problem. The only case that is reasonable to complete is the Toda 
system of type G 2 . In that case the second Poisson bracket should be a homogeneous 
bracket of degree 5 (a conjecture of Flaschka states that the degrees of the independent 
Poisson tensors coincide with the exponents of the corresponding Lie group). The other 
exceptional cases are even more complicated. It is a nontrivial task even to write down 
an explicit Lax pair for the systems and therefore the methods of this paper will be 
difficult to apply. 



• Full Kostant Toda 

One has a tri-Hamiltonian formulation of the A n system but no hierarchy. In the case 
of generalized full-Kostant Toda lattice (associated to simple Lie groups) one could 
seek to find similar structures as in the present paper. So far, nothing is known. The 
interesting feature of these systems is the presence of the rational integrals that are 
necessary to prove integrability. The Lie-algebraic background of the systems will 
certainly play a prominent role. 

• The Volterra or KM-system The multi-Hamiltonian structure of this system was 
first obtained in [64]. However, there is a symplectic realization of the system, due to 
Volterra, and it would be interesting to find a recursion operator in that symplectic 
space that projects onto the known hierarchy (as in section 5.1). 

• Bogoyavlensky- Volterra lattices 

There is also an interesting connection with the corresponding generalized Volterra 
systems also defined by Bogoyavlensky [65] in 1988. It seems that the multiple Hamil- 
tonian structures of the Volterra and Toda lattices are in one-to-one correspondence 
through a procedure of Moser. Multiple Hamiltonian structures for the generalized 
Volterra lattices, were constructed recently by Kouzaris [66], at least for the classical 
Lie algebras. The relation between the Volterra systems of type B n and C n and the 
corresponding Toda B n , C n systems was demonstrated in [67]. The connection between 
Volterra D n and Toda D n is still an open problem. 

• Independence of Poisson structures Equations (53) and (150) are remarkable; 
one Hamiltonian system, infinite formulations. On the other hand the systems are 
finite dimensional and after some point some dependencies should occur. Indeed, the 
integrals are not all independent. It is natural to ask a similar question about 
the infinite sequence of Poisson structures. Do they become dependent after a certain 
point? Unfortunately, there exists no widely accepted definition of independence for 
Poisson tensors. 

• Is the Toda lattice super— integrable? This conjecture should be true. A number of 
well-known systems are super-integrable, i.e. the free particle, the harmonic oscillator 
and the Calogero-Moser systems. In the case of the open Toda lattice, asymptotically 
the particles become free as time goes to infinity with asymptotic momenta being the 
eigenvalues of the Lax matrix. Therefore, the system behaves asympotically like a 
system of free particles which is super-integrable. 
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